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ABSTRACT 

We present the results of an archival XMM-Newton study of the bright X-ray point 
sources (L^ > 10 38 ergs“ 1 ) in 32 nearby galaxies. From our list of approximately 100 
point sources, we attempt to determine if there is a low-state counterpart to the Ultra- 
luminous X-ray (ULX) population, searching for a soft-hard state dichotomy similar to 
that known for Galactic X-ray binaries and testing the specific predictions of the IMBH 
hypothesis. To this end, we searched for “low-state” objects, which we defined as objects 
within our sample which had a spectrum well fit by a simple absorbed power law, and 
“high-state” objects, which we defined as objects better fit by a combined blackbody 
and a power law. Assuming that “low-state” objects accrete at approximately 10% of 
the Eddington luminosity (Done & Gierlinski 2003) and that “high-state” objects ac- 
crete near the Eddington luminosity we further divided our sample of sources into low 
and high state ULX sources. We classify 16 sources as low-state ULXs and 26 objects as 
high-state ULXs. As in Galactic black hole systems, the spectral indices, T, of the low- 
state objects, as well as the luminosities, tend to be lower than those of the high-state 
objects. The observed range of blackbody temperatures for the high state is 0.1-1 keV, 
with the most luminous systems tending toward the lowest temperatures. We therefore 
divide our high-state ULXs into candidate IMBHs (with blackbody temperatures of ap- 
proximately 0.1 keV) and candidate stellar mass BHs (with blackbody temperatures of 
approximately 1.0 keV). A subset of the candidate stellar mass BHs have spectra that 
are well-fit. by a Comptonization model, a property similar of Galactic BHs radiating in 
the “very-high” state near the Eddington limit. 
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Subject headings: galaxies: general — surveys - — X-rays:binaries — accretion, accretion 
discs 


1. Introduction 

Through X-ray observations of nearby galaxies, a class of Ultraluminous X-ray (ULX) sources 
has emerged. These are pointlike, non-nuclear sources with observed X-ray luminosities greater 
than 10 39 ergs _1 (Miller & Colbert 2004). Of most interest are those sources with bolometric 
luminosities in excess of the Eddington limit for a 20 M® black hole, or L boi > 2.8 x 10 39 erg s -1 . 
The true nature of these sources is unclear, and this class most likely includes several, different types 
of objects. Though some of these sources are located within a few parsecs of their host galaxy’s 
dynamical center, they do not exhibit many of the characteristics of active galactic nuclei (AGN). 
Because the ratio of X-ray to optical flux is a factor of 10 greater than that of AGN (Anderson et 
al. 2003; Stocke et al. 1983), these objects are fairly easy to recognize in X-ray imaging data. 

Assuming that the Eddington limit is obeyed by black hole accretion, the existence of such 
luminous non-AGN sources presents a puzzle. Several models have been proposed to account for 
the high luminosities of the ULXs. Among these are relativistic and non-relativistic beaming from 
stellar-mass black hole systems (Reynolds et al. 1997; King et al. 2001; Kording et al. 2002) and 
accretion of matter into intermediate mass black holes (IMBHs) . In several ULX systems (NGC 
1313 X-2, M81 X-9, etc.), detection of emission nebulae surrounding the ULX supports isotropic 
emission from the central source (Pakull & Mirioni 2003), which cannot be described through 
beaming. Further, a number of ULX (NGC1313 X-l, etc.) X-ray spectra are best fit with combined 
multi-component blackbody (MCD) and power law fits, similar to Galactic black holes in their high- 
state. Recently, Miller, Fabian, & Miller (2004) find that many spectral fits of ULXs require cool 
accretion disk temperatures of approximately 100 eV. The theoretical relationship between black 
hole mass and disk temperature (T oc M -1 / 4 ) has been observed to hold true for stellar mass 
(typically around IkeV) and supermassive (around 10-100 eV) black holes (Makishima et al. 2000; 
Porquet et al. 2004; Gierlihski & Done 2004). Using these scaling relations, the cool accretion disk 
ULXs would correspond to a population of high-state IMBHs with masses of % 16 — 10 4 M®. 

If ULXs do not obey the Eddington limit, they could be the result of an outburst (such as 
can occur in low mass X-ray binaries within our own Galaxy). Jonker & Nelemans (2004) find 
evidence for approximately 5 Galactic black hole X-ray binaries which exhibit luminosities in the 
ULX range during outbursts. These sources would appear as transient ULXs. The typical time 
scale for outburst of Galactic X-ray transients is a few days to rise from quiescent level with a 
decline from peak brightness to quiescent value of 30 - 40 days (Chen, Schrader, & Livio 1997). 
Another possible explanation is super-Eddington emission from accretion disks surrounding stellar 
mass black holes (Begelrnan 2002; Ebisawa et al. 2003). Sources of this type would be expected 
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to have soft X-ray components well modeled by hot accretion disks -(& 1.3 - 2.0 keV) similar to 
superluminal X-ray sources in the Galaxy (e.g. Belloni et al. 1997). 

Likely, ULXs include a variety of different objects with both isotropic and non-isotropic emit- 
ters. However, if some ULXs do indeed represent a class of high-state IMBHs, similar to the 
high/soft (thermal dominated) state stellar mass black holes in our galaxy, we might also expect to 
see the low-state objects from this same population. In Galactic black hole systems, the low-state 
is generally characterized by lower luminosity, with L< 0.1 L Edd. (Done & Gierlinski 2003), and a 
power law photon spectrum, typically with index T « 1.7 (McClintock & Remillard 2004). Indeed, 
the existence of some ULX sources (IC 342 X-l, NGC 5204 X-l) as possible low-hard (pure power 
law) state IMBHs, well-fit by simple absorbed power laws, have been noted from Chandra obser- 
vations by Roberts et al. (2004). In this study we seek to test a direct prediction of the IMBH 
hypothesis; namely, whether there is a class of sources with properties consistent with what we ex- 
pect of low-state IMBHs. This requires two major assumptions: (1) that the emission from ULXs is 
isotropic and (2) that IMBHs exhibit the same states (whose classification was based on luminosity 
and spectral form) as stellar mass black holes. Our goal is to find these “low-state” sources, if they 
exist, classify the properties of both high-state and low-state ULXs, and test whether these data 
are consistent or inconsistent with the predictions of the IMBH hypothesis. 

We present the results of a detailed analysis of ULXs in nearby galaxies observed with the 
European Space Agency’s XMM-Newton observatory. Only XMM-Newton provides the count rates 
and bandpass necessary to distinguish different spectral models for most ULXs, accurately deter- 
mine both the temperature of the thermal component expected for high-state objects, and de- 
termine whether this component is required in the spectral modeling of these objects. Since the 
XMM-Newton X-ray spectra of ULXs are similar in quality to spectra for Galactic X-ray binaries 
obtained in the 1980s, our spectral classification in this paper will remain purely schematic. Thus, 
our classifications as low and high state objects are a first approximation, based on the quality of 
the spectra available. 

In Section 2, we detail the observations examined from the XMM-Newton archives and explain 
the data analysis for the individual point sources. In Section 3, we discuss the spectral fitting tech- 
nique as well as simulations we conducted to determine their validity. We discuss the implications 
of our results in Section 4. 


2. Observations and Data Reduction 

The data used in this investigation were drawn from the XMM-Newton public data archive. 
Assuming that low-state ULXs exist in the luminosity range of It) 38-39 ergs -1 , we conducted sim- 
ulations to determine the optimum criteria for observations capable of resolving point sources of 
this luminosity. Our simulations provided a guide for choosing which of the vast number of archival 
XMM datasets we should examine. This luminosity range was chosen on the assumption that an 
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approximately 100 black hole would radiate at « 10% of the L Edd in the low-state (Done & 
Gierlinski 2003). 

Within the luminosity range of interest (Ljy > 10 38 erg s -1 ), there are a number of known 
objects that could be confused with ULX sources. One type of source is supernova remnants (SNRs). 
These sources are often easy to distinguish based on their characteristic spectrum: with poor signal 
to noise we expect a steep power law and as the signal to noise increases, emission lines become 
clearly visible. Super-Eddington accreting neutron stars have been observed to have luminosities 
within this range for a short period of time. Neutron star X-ray binaries often have spectra well 
fit by a hot multi-color disk blackbody model, or with low signal to noise, by a bremsstrahlung 
model. Both models have similar curvature and a 0.7-2. 0 keV blackbody model is indistinguishable 
from the bremsstrahlung model. We chose to use the bremsstrahlung model because it is the 
simpler model and gives an adequate qualitative description to the data. We expect that for low 
temperature bremsstrahlung sources, the spectrum should be easily distinguishable from a power 
law with T 1.7 (as is expected for a low-state object). If, however, the neutron star spectrum 
has kT > 5keV, as observed for some NS X-ray binaries, our simulations show that we can not 
distinguish between the power law and bremsstrahlung models. 

The most common sources we expected to find in this luminosity range were the analogs to 
Galactic black hole X-ray binaries in their high-soft (thermal dominated) state. These sources 
typically have spectra well fit by a blackbody with temperature of ps 1.0 keV combined with a 
power law with index T ps 2.5. Our simulations sought to determine the number of photons 
required to distinguish between spectral fits corresponding to a power law model with T ps 1.7 and 
a combined blackbody and power law model. These models qualitatively correspond to those of a 
low-state (pure power law spectrum) and high-state (thermal dominated spectrum) X-ray binary. 
Since we do not know the proper normalization between the blackbody and power law components 
for high-state objects (it varies from source to source), we tested whether each of the components 
separately, e.g. blackbody or a steep power law, could be distinguished from the simulated “low- 
state” spectrum. We chose to simulate spectra in XSPEC using the command fakeit none. We used 
generic response and ancillary response matrices. Simulating a power law model with a T = 1.7, we 
found that for 200, 400, and 1000 counts, these models were distinguishable at > 99% confidence 
from a blackbody source (with kT constrained to the range of 0.6 to 1.3 keV, similar to that of 
Galactic black holes). We found that for a lower number of counts the distribution in T values 
increases to include a larger range of T values (i.e. T = 1.3 - 2.0 compared to T — 1.5 - 1.7). 
Simulating a power law with T = 2.5, we find the same trend. We determined that at roughly 400 
counts the distributions of T from a T — 1,7 and T = 2.5 power law become entirely separable at 
> 99 % confidence. 

In order to distinguish between the different spectral fits for objects with Lx ~ 2 x 10 38 ergs -1 , 
we select all galaxies that were observed for at least 10 ks (with the exception of the bright ULX 
in NGC 5408, which had enough photons for analysis despite the low exposure time) with XMM- 
Newton and that are closer than 8 Mpc. We estimated that these criteria would give us a minimum 
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of 400 counts for objects with Lx > 2 x 10 38 ergs -1 . We emphasize that the criteria quoted, based 
on the simulations, were used as a guide in choosing the sample of galaxies examined in this study. 
These simulations are not used as the statistical basis for our object-by-object analysis (discussed 
in Section 3). 

Our sample of galaxies is selective in that it represents objects of interest in the X-ray band. We 
include details on these host galaxies in Table 1 . XMM- Newton spectral information of individual 
X-ray sources had previously been published for approximately 60% of the host galaxies. We include 
references in the alternate ID column and footnotes of Table C6. We do not compare our results 
with these previous studies on a source by source basis. 

We found that abstracts describing the proposals for XMM-Newton observations were available 
for only 13 of the 32 galaxies examined. Of these 13, only one observation cited the motive as a 
study of ULXs (NGC 1313). However, 7 of the remaining 19 galaxies contained sources classified 
as IXOs, intermediate X-ray objects, by Colbert & Ptak (2002). If the remaining galaxies were not 
studied due to their ULX population, the effects of bias are small with roughly 25% of the sources 
studied explicitly due to their connection with ULX sources. Our host galaxies include only spirals 
and irregulars. Figure 1 displays the distribution of galaxy type. 

We reduced the data using the XMM-Newton Science Analysis System (SAS) version 6.0.0. 
Since the processed pipeline products (PPS) were created with earlier versions of SAS, the observa- 
tion data files (ODF) were used to produce calibrated photon event files for the EPIC-MOS and PN 
cameras using the commands emchain and epchain. Following this, the events tables were filtered 
using the standard criteria outlined in the XMM ABC Guide. For the MOS data (both MOSl and 
MOS2 cameras), good events constitute those with a pulse height in the range of 0.2 to 12 keV and 
event patterns that are characterized as 0-12 (single, double, triple, and quadruple pixel events). 
For the PN camera, only patterns of 0-4 (single and double pixel events) are kept, with the energy 
range for the pulse height set between 0.2 and 15 keV. Bad pixels and events too close to the edges 
of the CCD chips were rejected using the stringent selection expression “FLAG == 0” . 

Time filtering was applied as needed by editing the light curve produced in xmmselect for 
the entire observation. Flare events (distinguished by their high count rate) detected in all three 
cameras, were cut using the tabgtigen task as outlined in the ABC Guide. Typical count rate 
parameters for filtering were ’RATE < 5’ for MOS detectors and ’RATE < 20’ for the PN detector. 
Such filtering was only done as needed. Pre-filtered exposure times are listed in Table 1. The 
number of counts from the filtered net exposure times for the individual sources are listed in 
Table C6. We note that the filtered data are not always sufficiently clean that a high signal-to- 
noise is maintained up to 10 keV. Sources with a high background flux level, relative to the source 
spectrum, show poorer signal-to-noise in the spectrum above lkeV. 

Before extracting spectra of the brightest sources, contour maps of the X-ray observation 
were overlaid on Digital Sky Survey (DSS) images. This ensured that bright foreground stars and 
background AGN were easily distinguished, and thereby not included in our sample. Also, we 



- 6 


checked the XMM-Newton source positions with NED and SIMBAD to determine if they coincide 
with any known background galaxies or QSOs. A list of these bright fore-ground or background 
sources is included in Table C7. 


3. Spectral Fitting 

Spectra for the bright point sources were extracted based on their apparent brightness in the 
CCD images. No explicit source detection algorithm was necessary. We used the SAS task especget. 
With this task we created spectra (for both the source and background), response matrices, and 
ancillary response files for all three EPIC cameras, when possible. The typical extraction radius 
was 20 arcseconds, but depending on both the size and proximity of a source to another source, the 
extraction radius ranged from 9 - 87 arcseconds. Background spectra were extracted either in an 
annulus centered on the source, or in a circle of appropriate size away from the source, depending 
on the proximity of the candidate source to other X-ray sources. Annuli were used for sources that 
were not located within a few arcseconds of another source, thus annular background extraction 
radii were not used for sources with small extraction radii. For sources in crowded regions, we used 
circular extraction radii close to the source. We extracted backgrounds close to the source in order 
to correct for emission local to the ULX. Once the spectra were obtained, they were rebinned to 
require at least 20 counts per bin, using the command grppha in LHEASOFT. The list of sources, 
with position and count information, is included in Table C6, We only included sources for our 
spectral studies that had 400 or greater PN counts (or MOS for sources in NGC 253 and M81 and 
NGC 4945 XMM3 and NGC 2403 XMM1, for which PN spectra were not available).. 

The extracted spectra were fit with standard models in XSPEC vll.3.1. For each source, we fit 
the PN and MOS spectra simultaneously in the 0.3-10 keV range. We allowed a free normalization 
constant to account for the differences in flux calibration between the three cameras (similar to 
Jenkins et al. (2004)). Each source was first fit with an absorbed single component model. In all 
cases we used the standard absorption model wabs , leaving the column density as a free parameter. 
For those sources where the hydrogen column density was unconstrained, we fixed the value to the 
Galactic foreground value listed in Table 1. 


3.1. Single- Component Sources 

Results of the single-component fits are seen in Table 2. We include in this table only the 
best-fit parameters for those sources best described by a single-component model. The addition of 
a blackbody component to these single-component fits changes the y 2 value by a negligible amount 
and therefore is not statistically significant. More specifically, the addition of a blackbody model 
to the power law fit corresponds to a Ay 2 < 2.3, which is the « 68% confidence level using the 
F-test for two degrees of freedom. 
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The flux values quoted represent the unabsorbed flux in the PN spectra, in the 0.3-10 keV 
band. All errors quoted, here and subsequently, correspond to the 90% confidence level for one 
degree of freedom (Ax 2 = 2.71) . The luminosities were calculated from the unabsorbed flux using 
the distances quoted in' Table 1. Both flux and luminosity correspond to those of the best fit model 
(power law or bremsstrahlung). It should be noted that since our selection criteria was based on 
a count rate cutoff, due to the variety of spectral forms, the inferred luminosity cutoff will not be 
uniform. 

In Table 2 we denote the single component model we choose as the better fit in bold. This 
notation also indicates the model (power law or bremsstrahlung) used to compute the quoted 
flux. For « 46% of the power law sources, the % 2 difference (< 2) between the power law and 
bremsstrahlung models is only marginally different. Of these sources, the average kT value for 
the bremsstrahlung fit is 5.54 keV. From our simulations we find that at high temperatures the 
bremsstrahlung fit becomes indistinguishable from a simple power law. Thus, given the high tem- 
peratures of the bremsstrahlung fits for these sources, they are equally well described by the power 
law model. Typical kT bremsstrahlung values for accreting neutron stars are from 3.0 to 7.0 keV 
(Jones 1977). 


3.1.1. “Low State” ULXs 

From these single-component sources, we classify 16 sources as low-state ULX sources. This 
classification is based on (1) the shape of the spectrum, well-fit by an absorbed power law, (2) 
the luminosity of the sources (they needed to be luminous enough to be included in our sample, 
Lx > 10 38 ergs _1 ), and (3) the X-ray location of the object within the optical galaxy (based on 
Digital Sky Survey (DSS) images) . The third criterion was important in limiting the effects of 
contamination from fore-ground and background sources within our sample of ULX sources. Thus, 
we overlaid X-ray contours from the XMM image on the DSS images, determining the location of 
the X-ray source as within or outside of the optical extent of the galaxy. We note that for two 
of the sources classified as low-state ULXs (Holmberg I XMM2 and NGC 2403 XMM4), there is 
uncertainty of whether the X-ray source is in fact within the optical galaxy due to the quality 
of the optical DSS images. The images used to determine the third criterion are available online 
(http://www.astro.umd.edu/~lwinter/second3.html). We further discuss these criteria in Section 
4. 

We state that the first criterion for classification as a low/hard state ULX is a spectrum well-fit 
by a power law. Of the 30 sources in Table 2, three sources were clearly not well-fit by either the 
power law or bremsstrahlung model. Seven of the remaining 27 sources were clearly not within the 
optical extent of the host galaxy. Of the remaining sources excluded from classification as a low- 
hard state object IC 2574 XMMl was better fit with a bremsstrahlung model (with A% 2 = 13.4). 
In an additional observation of the source NGC 4258 XMM2, a simple power law model is not an 
adequate fit to the data (while the luminosity of the source is not within the “ULX” classification 
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range mentioned in the introduction) . For the remaining sources, NGC 247 XMM2 and NGC 253 
XMM2, there was sufficient doubt on the spectral form where the bremsstrahlung and power law 
model as well as the addition of a thermal component all yielded adequate fits to the spectra. 
Therefore, we excluded these sources from a classification, noting the ambiguity of the model fits 
for these sources. 

For those sources we classify as ULXs, we include computed bolometric luminosities in Table 5. 
To compute the bolometric luminosities for these ULX sources, we used the exponentially cutoff 
power law spectrum of Magdziarz & Zdziarski (1995) , model cutoffpl in XSPEC , with a cutoff 
energy of lOkeV. From observational studies of Galactic X-ray binaries, it has been observed that 
low-state objects have spectra that cut-off at high energies (> 10 — 200 keV) (Zdziarski & Gierlinski 
2004). Thus we chose the exponential model cutoffpl over a simple power law. This also minimizes 
the total luminosity for fiat power law sources. We computed an unabsorbed bolometric flux in 
the 0.1 - 100 keV range through use of the dummyresp command (which extends the model beyond 
the observation’s energy range) . The luminosity was then computed using the distances listed in 
Table 1. We quote these values as L cutofjpl (the luminosity obtained from extrapolating the power 
law portion of the spectrum as an exponentially cut-off power law) in Table 5. We note that these 
values represent an upper limit on the bolometric luminosity for steep power law (F > 2) objects, 
since we would expect the power law component to cutoff at some low energy. However, for shallow 
spectrum (P < 2) sources L cutoffpl is a lower limit. This is because, schematically, a steep power 
law diverges at low energies while a shallow power law diverges at high energies. 


3.2. Two- Component Sources 

For a number of sources, we found that an improvement in reduced x 2 was achieved through 
fitting their spectra with an absorbed two-component blackbody and power law model. We chose a 
simple blackbody model over the multi-component disk model, diskbb, for purely schematic reasons. 
Namely, observations of galactic X-ray binary systems were fit with blackbody models in the 1980s, 
when the signal-to-noise of these objects was comparable to that for our XMM data for ULX sources. 
We also note that the dislcbb model does not give an entirely accurate physical description of the 
data as it neglects the effects of general relativity. As a schematic model, the blackbody model is 
simpler than diskbb, with the same number of degrees of freedom. In addition, for low temperatures 
both models yield virtually identical temperatures. For this study, we chose the simpler model. We 
defer to a further paper a discussion of the different models for the thermal component. 

In Table 3 we present the results for the sources which are fit significantly better by the two- 
component model, these are sources where the improvement in x 2 is greater than 8 (determined from 
our simulations in Appendix A). We include in Table 3 the improvement in x 2 of the two-component - 
fit over the simple power law. We include the power law best fits to these sources in the appendix 
for comparison with other analyses. In order to determine whether the blackbody component is 
statistically significant for the sources fit with a two-component model, we simulated spectra based 
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on accurate modeling of some of the brightest sources: NGC 247 XMM1, NGG 5408 XMM1, and 
Holmberg II XMM1. These sources span the observed range of the ratios of the blackbody to power 
law component and thus represent those from our sample with a weak blackbody relative to the 
power law component, intermediate case, and a strong blackbody, respectively. Our simulations are 
described in full in Appendix A. We found that, using a A^ 2 > 8 criterion, which corresponds to the 
99% significance level as according to the F-test for the addition of two extra parameters, we can 
readily detect the strong and intermediate thermal components in all spectra with more than 400 
counts. The weak thermal emission cannot be detected in 400 count spectra, but is readily detected 
in 2000 count spectra. This gives us confidence that our results are statistically meaningful. 


3.2.1. “High-State ULXs” 

Of the sources in Table 3, we classified high-state (or thermal dominated) ULX sources based 
on three criteria: (1) spectra characterized by an absorbed power law and blackbody model, (2) 
luminosity, and (3) X-ray source within the optical extent of the host galaxy. The luminosity 
criteria required that these sources have unabsorbed luminosities Lx > 3 x 10 39 ergs -1 (we used 
Lx = 2.7 x 10 39 ergs _1 as our hard cutoff). If the sources are radiating at the Eddington luminosity, 
this cutoff luminosity corresponds to objects with masses greater than 20 M®. 

From Table 3, 27 observations are recorded with Lx > 2.7 x 10 39 ergs -1 . The addition of a 
thermal component to these sources is statistically significant over a pure power law model. Of the 
27 observations, 3 correspond to multiple observations of a single source. From an analysis of the 
DSS images, all 24 of these sources are within the optical extent of their' host galaxies. However, 
M51 XMM5 appears to be coincident with the center of its host, a dwarf companion galaxy to 
M51. The location, coupled with the high luminosity (Lx — 1-9 x 10 42 ergs -1 ) leads us to classify 
this source as an AGN. We also excluded two sources (NGC 1313 XMM2 and M81 XMM2) from 
our sample of high/soft state ULXs due to their previous identification as supernovae. Of the 
remaining 21 sources, NGC 253 XMM4 had a luminosity of 2.5 X 10 4 ° ergs -1 in one observation 
and 2.2 x 10 39 ergs -1 in a second. This significant change in luminosity, with one observation below 
our luminosity cutoff and another a factor of 10 higher than the other, led us to exclude this 
source as a high/soft state ULX. It is likely that this source is a stellar mass X-ray binary within 
its host galaxy, where one of the observations captured the source in an outburst. 

In Table 4, we list sources that have Ax 2 values less than 8 for a single observation. Most 
of these sources have weak blackbody normalizations compared to the power law normalization. 
We classify these sources as being well-fit by a two-component model while acknowledging the 
uncertainty in the fit as determined by the simulations. The addition of the thermal component 
is not significant enough for these sources to be classified with certainty in either Table 2 or 3. 
The simple power law fits for these sources are included with those for sources in Table 3 in the 
appendix. We note that due to their high luminosity we included six of these sources (NGC 4490 
XMM2, NGC 4490 XMM3, NGC 4736 XMM1, M51 XMM2, M51 XMM6, and M101 XMM3) with 
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uncertain fit parameters as ULX high-state sources. Two of these sources (NGC 4490 XMM3 and 
M51 XMM6) had unabsorbed luminosities > 10 times the 3 x 10 39 ergs -1 cutoff used for high-state 
ULX classification. The other four sources had luminosities above the threshold, as well as weak 
blackbody components compared to the power law (see the appendix for simulations). We used 
these points to justify including these sources with the Table 3 sources in the following discussions 
with the proviso that their spectral fits do not indicate absolutely the necessity of the additional 
thermal component. For this reason, we denote these sources with a special symbol (a circle) in 
subsequent figures while including them as “high-state” ULX objects. 

For our ULX sources modeled by a combined blackbody and exponentially cutoff power law, we 
computed bolometric luminosities using two methods. The first method is recorded as L cuto f f p i in 
Table 5. We computed the flux from 0.1 - 100 keV using an unabsorbed blackbody and exponentially 
cutoff power law model using the XSPEC command dummyresp. For the second method, recorded 
as L bol in Table 5, we estimate a more accurate bolometric luminosity calculated from the flux 
in the range of 2x kT- 100 keV where kT is the blackbody temperature obtained from the model. 
In Galactic X-ray binary systems, the power law component of the X-ray spectrum is believed 
to be from Comptonization in a corona. The photons supplying this energy originate from the 
blackbody continuum emanating from the accretion disk. Thus, a natural cutoff for this power 
law component occurs at the peak emission of the blackbody (which is approximately 3x kT). The 
estimated values (obtained from cutting off the combined unabsorbed blackbody and cutoff power 
law model at the value 2x kT) differ with regard to the full estimate (flux from the fully integrated 
blackbody added to the separate flux from the cutoff power law from 3x kT to 100 keV) depending 
on the normalization factors used (for both the blackbody temperature and the spectral index T). 
Choosing three sources displaying a range of blackbody to power law strength (Holmberg II XMM1, 
NGC 253 XMM1, and IC 0342 XMM3) we found that the estimated values were within 88.3, 95.1, 
and 96.8% of the more complete estimation. Given their dose proximity (within approximately 
90%) we quote these estimated values as a good approxmiation to the bolometric luminosity. 

We note that our bolometric luminosities for all of the classified ULX sources, on average, are 
a factor of 1.08 greater than the X-ray luminosities in the 0.3 — lOkeV band for the objects best 
fit by a combined blackbody and power law. Thus, to good approximation, the X-ray luminosity 
is the bolometric luminosity. However, for the objects best fit by a simple power law (low-state 
ULX sources) , the average bolometric luminosity is roughly a factor of 7 greater than the X-ray 
luminosity in our band. This average is dominated by the steep power law objects, in particular 
Holmberg II XMM1 (F — 3.09). Excluding this object, we get an average bolometric luminosity 
that is 2.8 times the X-ray flux and more indicative of the general properties of these power law-fit 
objects. 



3.3. Additional Sources 


In addition, in this large sample of point sources, we came across a number of objects whose 
spectra were not well fit by the models we employed. These sources have luminosities exceeding 
Lx ~ 10 38 erg s -1 , if they are associated with the host galaxy, and are placed in Tables 2 and 3 
as well as Appendix B. These sources include two supersoft sources, one possible AGN, and three 
sources well fit with additional absorption models (including a partial covering model and a model 
of hot gas) . We briefly describe these sources in the appendix. 


4. Discussion 

We have determined best-fit spectral parameters of the bright X-ray sources in 32 nearby 
galaxies. In choosing three “standard” models for our study, we hoped to accurately separate 
high and low state ULXs from other types of luminous X-ray sources. We specifically chose to fit 
the data with the bremsstrahlung model in order to identify neutron star X-ray binaries within 
our sample. The models we used are purely schematic, and they do not physically explain the 
phenomena occurring, but are standard and qualitatively simple models often used to fit the spectra 
of Galactic X-ray binaries. 

We cross-referenced the X-ray positions of our sources with both NED and SIMBAD in order to 
identify known supernovae, galaxies, and stars. In addition, we examined the DSS optical images to 
place the position of our sources within their respective galaxies. Such analysis aimed to minimize 
contamination of our sample of ULXs with bright background and foreground sources. 

Further, we examined XMM-Newton’s Optical Monitor data in the visual bands (U, B, V). 
The XMM PPS contain point source detection files for the OM data. We overlaid these point 
source detections with X-ray contour maps in order to determine the brightest possible optical 
count rates for the X-ray sources, which were then converted into fluxes using the OM calibration 
documentation. In Figure 2, we plot the distribution of the logarithm of the X-ray to optical flux 
for the brightest possible optical counterpart inside the XMM- Newton error circle. Only 13 of the 
32 host galaxies had visible band OM data during the observations. Of these 13 galaxies, 40 of the 
X-ray sources were in the range of the OM data and only 14 were coincident with an optical point 
source. Therefore, the majority of our sources have X-ray/optical flux ratios that are larger than 
those displayed. Figure 2 illustrates the lowest possible X-ray/optical flux ratios and also, by the 
sparsity of sources included in the diagram, it illustrates the fact that a majority of the sources have 
no obvious optical counterpart and thus have very large X-ray/optical flux ratios. We estimate the 
point source detection limit of the OM U filter as approximately 1.24 x 10 -14 erg cm -2 s -1 . For an 
unabsorbed X-ray flux of 1.0 X 10~ 12 ergcm -2 s -1 , typical of objects with Lx ~ 2 x 10 38 ergs -1 
located at a distance of 8Mpc, this corresponds to log {f x / fopt) ■— 1-9. Therefore, the average value 
for our sources should fall around 2 or greater. The average distribution for QSOs and AGN centers 
around 0 and 0.8 for BL Lacs (Anderson et al. 2003). Our objects have ratios of L x /L opt at least 
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10 times higher than those of AGN and 100 times greater than stars. 

Recently, Gutierrez & Lopez-Corredoira (2005) identify six ULXs from the catalog of Colbert 
& Ptak (2002) as QSOs. They hypothesize that a large number of ULXs may in fact be quasars at 
higher redshift than their supposed host galaxy. However, unlike the objects studied in Gutierrez & 
Lopez-Corredoira (2005) , our ULX sources are all spatially coincident with the optical host galaxy. 
In addition, a majority of our ULXs are not in the proximity of a noticeable optical point source. 
The X-ray/optical flux ratios of our sources are much larger, on average, than might be expected for 
a QSO. It is also worth noting that while some cataloged ULXs may be QSOs, optical identifications 
have been made associating other ULXs with a type B supergiant companion (Kuntz et al. 2005; 
Liu et al. 2004). 


4.1. Classification Criteria 

The spectral fits indicate that to high statistical probability (see appendix A) we can distinguish 
a class of low-state ULXs from the high-state objects. This is assuming, as indicated in the 
introduction, that ULXs are isotropic emitters with luminosity and spectral form similar to Galactic 
stellar-mass X-ray binaries. In section 3, we stated that our ULX classification depends upon three 
criteria: (1) spectral form, (2) luminosity, and (3) location of the X-ray source within the optical 
host galaxy (as determined from the DSS images). We have chosen simple, parametric ‘non- 
physical’ models for the spectra because the signal to noise of most of the observations does not 
allow anything else to be constrained. 

Of the sources in Table 2, 16 are classified as “low-state” objects or low/hard state ULXs, 
having unabsorbed luminosities > 10 38 ergs -1 and spectra that are well fit by power law models. 
Throughout this paper, we use the term low-state ULXs to include “low-state IMBH candidates” 
(sources with Lj < 3 x 10 39 erg s'" 1 and Spectra well-fit by a simple absorbed power law) and low- 
state sources with luminosities that clearly classify them as ULX sources ( Lx > 3 X 10 39 ergs -1 ). 
These low-state ULX sources are listed in Table 5. 

In the Spectral Fitting section, we noted that a power law and high temperature bremsstrahlung 
model are indistinguishable. Therefore, it is important to consider the luminosity of these sources 
in the claim that they are not neutron star X-ray binaries accreting at the Eddington luminosity. 
Of the low-state ULX sources, only two of the 16 sources have bolometric luminosities below the 
Eddington luminosity of a 3 Mq object (« 4 X 10 38 ergs -1 ), corresponding to the maximum mass 
of a neutron star. All of the sources have values exceeding the Eddington limit for a 2 Mq neutron 
star. 

Further, 26 sources have unabsorbed Lx > 3 x 10 39 ergs -1 , corresponding to L « Lem at M 
> 20 Mq as expected for “high-state” IMBHs, and spectra that are well fit by combined blackbody 
and power law models. These are “high-state” objects. The spectral fits for these sources are listed 
in Tables 3 and 4. In a statistical sense, we find that the greater the number of counts in the 
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observation the greater our confidence in the thermal component contributing to a better fitting 
model. We explain our confidence levels obtained from spectral simulations in the appendix. 

In addition to these high and low state ULXs, we find a large number of sources best fit by 
a combined blackbody and power law model but below our threshold of Lx « 3 X 10 39 ergs -1 for 
a high-state ULX (listed in Table 3). Many of these sources may be accreting stellar mass black 
holes with M< 20 Mq. Some of these non- “ULX” sources were found away from the optical extent 
of the targeted galaxy (from our analysis of the DSS images), and therefore may be background 
AGN. 


4.2. Low-State ULXs 

For Galactic black hole X-ray binaries, spectral indices of low-state (or power law dominated) 
objects are typically lower than those of high-state objects, with < F >i ow & 1.7 and < F >hi g h~ 2.5 
(McClintock & Remillard 2004). In Figure 3, we plot the distribution of the spectral index for both 
high-state and low-state objects. The spectral index for the high-state objects is the value of T 
from the two-component fit. As in the Galactic sources, it is clear that the spectral indices of the 
high-state objects are indeed larger. Of further interest, the distribution of spectral index for low- 
state objects looks remarkably similar to the distribution of spectral index for moderate luminosity 
quasars, many of which are thought to be the analogs of low-state black holes (Porquet et al. 2004). 
This supports the classification of these objects as accreting black holes. 

For the high-state objects, we find mean values of F = 2.46, with a root mean square (rms) 
deviation of S= 0.12, and Lx = 1.4 X 10 4O ergs -1 , log(S) = 1.6. This calculation excludes the 3 
objects with spectral indices greater than 3.5. For the low-state objects, we find mean values of 
F = 2.09, with a rms deviation of S= 0.10, and Lx = 2.2 x 10 39 ergs _1 , log(5) = 2.1. This value 
of F = 2.1 is softer than the typical hard-state value of « 1.7, but within the 1.5 < T < 2.1 range 
used to classify this state for Galactic X-ray binaries (McClintock & Remillard 2004). Computing 
a Kolmogorov-Smirnov two-sample test, separating the sources into the category of low-state or 
high-state, we find a likelihood of approximately 0.03 that the spectral indices belong to the same 
distribution. 

The low hard X-ray state of X-ray binaries is associated with a low accretion rate from the 
companion object with L< 0.1 L Edd (Done & Gierlinski 2003). Therefore, on average, we expect 
the luminosities of the low-state objects to be lower than the high-state objects. Figure 4 displays 
the luminosity of the objects as a function of the spectral index. On average, the highest luminosity 
low-state objects have luminosities lower than those of the high-state objects. 

The lower Lx values of the low-state objects imply that they may indeed be accreting at a 
lower rate than the high-state objects. This can further be seen in the bolometric luminosities 
listed in Table 5. If these objects are accreting at a rate similar to Galactic low/hard state black 
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holes (0.1 x Lgdd) (Done & Gierlinski 2003), we can estimate their masses as 

M _ Lm 
M & 0.1 x Lsdd 

with Leu as the Eddington luminosity for a 1 M® object (1.3 X 10 38 ergs _1 ). Our mass estimations, 
based upon our limits to the bolometric luminosities, yield masses of 20 — 1500M® (see Table 5), 
precisely what we might expect for a population of IMBHs. 


4.3. High- State ULXs 

If the high-state (thermal dominated) ULXs represent a class of intermediate mass black hole 
systems, their X-ray spectra should be well described by a combined blackbody and power law 
model. Scaling for the mass of the black hole, we would expect a relationship of Toe M -1 / 4 between 
black hole mass and blackbody temperature (Makishimaet al. 2000). This would indicate a thermal 
component of ~ 100 eV for masses of ~ 10 3 M®. A few objects have been reported to display this 
property (Miller et al. 2003; Roberts et al. 2005). In Figure 5, we graph the distribution of the 
thermal component for our classified high-state objects. 

We find that there are two peaks in the distribution among the thermal component, one at 
approximately 100 eV and another centered close to 1 keV. This could indicate two different classes 
among the high-state objects. It is possible that those objects with blackbody components near 
100 eV are indeed high-state intermediate mass black holes. We note that the soft excess in PG 
Quasars has also been modeled as a blackbody with kT S0 ^ ss 100 eV, but it has been suggested 
that this could be the result of a process not directly related to black hole accretion (such as the 
presence of a warm absorber: Gierlinski & Done (2004)). Another possible explanation is that the 
soft component is the result of ionized reflection from the disk (Ross & Fabian 2005). While the 
possibility exists that the “thermal” component of these 100 eV sources is not directly related to 
black hole accretion or is related in a “non-thermal” (i.e. ionized reflection) sense, as may be the 
case with the soft excess in PG Quasars, we assume that the soft component for the objects we 
classify as high-state ULXs originates from a thermal disk. We use this assumption to test the 
IMBH hypothesis, thus speculation on the nature of the soft component is beyond the scope of this 
paper. 

The second peak, centered around 1 keV, has a temperature reminiscent of the Galactic black 
hole systems in our own galaxy. These systems may thus be stellar-mass black holes accreting 
matter near the Eddington limit. If this were the case, we would expect the luminosities of the 
sources exhibiting a higher blackbody temperature to be lower than those with cooler blackbody 
components. In the second graph of Figure 5, we plot the relationship between blackbody tem- 
perature and Lx in the 0.3 - 10 keV band. Once again, two groups are seen in the distribution 
of high-state ULXs. The most luminous objects are those with low blackbody temperatures. On 
average, the less luminous sources exhibit higher blackbody temperatures. For the sources with 
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Lx > 1 0 40 erg s *, the mean blackbody temperature is 0.31 while the sources below this luminosity 
threshold have a mean blackbody temperature of 0.61. 

The second, low-luminosity, group in the distribution of high-state ULXs is clearly distinguish- 
able in both plots of Figure 5. We found that, with the exceptions of NGC 253 XMM1, M81 XMM1, 
and NGC 5204 XMMl, the spectra of these objects could also be well-described by an absorbed 
Comptonization ( compST) model (Sunyaev & Titarchuk 1980) used to fit galactic black holes in the 
“very high” state when they are radiating at the Eddington limit. This model simulates Compton 
scattering of cool photons on the hot electrons of a completely ionized plasma. We present the 
best-fit parameters for the Comptonization model in Table C9. 

This “very high” state has been observed (Miyamoto et al. 1991) in a few Galactic black holes. 
Yet another rubric for the very high state emerged in Kubota et al. (2001) and Kubota & Makishima 
(2004), where they identify this as the “anomalous” state, a state whose spectrum can be well fit 
by a Comptonized scattering model. Regardless of the name, our best-fit Comptonization sources 
likely fit into this category. The luminosities of these sources suggest that they are stellar mass 
black hole systems in this anomalous/ very high state. 

As with the low-state, we include mass estimates for our high-state objects in Table 5. We 
assume that the high-state objects are radiating at L Edd, resulting in a minimum mass if there is no 
beaming. We find masses of 1.6 — 38 M®, consistent with “normal” stellar mass BHs, for the sources 
well fit by the Comptonization model. The other high-state ULXs masses range from 17 - 1350 M® 
based on Eddington rates, analogous to the low-state ULX masses computed. 

It is important to note that the initial simulations (appendix) and discussions in the Spectral 
Simulations section need to be considered in relation to the impact they pose to our classification 
scheme and the results presented in these Sections. While it is indeed possible that some of the 
objects with a weak blackbody component and a relatively small number of counts would be mis- 
categorized as a pure power law spectrum, one can ask what such a possible situation would do 
to the correlations that we have seen. These putative objects, by assumption would have lower 
luminosities, however their temperatures are unknown and it is entirely unclear if they would 
destroy the kT/L(x) correlation. As we have shown in our simulations it is unlikely that the fitted 
power law index would change and thus the presence of a low state as indicated by the spectral 
index would not change. This would create a new type of object, one with a flat power law and 
a black body component, which is not seen in the Milky Way, nor among the high signal to noise 
objects. 


4.4. Temperature Gap 

In addition to the existence of ULXs with low blackbody temperatures, the temperature dis- 
tribution of the high-state (thermal dominated) ULXs (Figure 5, left panel) displays a “gap” which 
is of particular interest — there is a complete absence of objects with temperatures in the range 
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0.26 keV to 0.50 keV. It is tempting to take this as evidence for a gap in the mass distribution of 
these accreting black holes. Since, for a given luminosity, we expect the temperature to vary as 
T <x I 1/4 M -1 / 2 , this factor of two gap in the temperature distribution translates into a factor of 
four gap in the black hole mass distribution. 

If this result is borne out by further study, it provides an important clue to the origin and 
evolution of intermediate mass black holes. One popular idea is that intermediate mass black holes 
formed from the collapse of massive Population III stars (Madau &, Rees 2001). Models suggest 
that Pop III stars with zero age main sequence (ZAMS) masses in the range 25-140M© and above 
260M© collapse to produce black holes (Heger & Woosley 2002) whereas in the range of ZAMS 
masses 140-260M©, pair-instability supernovae lead to the complete disruption of the stars (i.e., 
no remnant black hole remains). Hence, this model for IMBH formation predicts a gap in the 
IMBH initial mass function in the range of approximately 60-200 M© (although this is uncertain 
on the low end due to the effect of the pulsational pair-instability on the pre-collapse core). One 
possibility is that the gap in our observed temperature distribution (and hence the inferred gap 
in the mass function) is due to this effect of the pair instability supernovae in Pop III stars. This 
would require that the current IMBH mass function is approximately the same as the initial IMBH 
mass function. In other words, it requires that most IMBHs (especially those just below the gap) 
have not grown significantly due to accretion since their formation and, hence, that the ULX phase 
represents a short fraction of the life-time of an IMBH (/ << f sa i/fH, where t sa . i « 45eo.i Myr is the 
e-folding timescale for Eddington limited black hole growth with radiative efficiency e = O.leo.i)- 

An alternative interpretation of the inferred mass gap is to suppose that two fundamentally 
different modes of formation lead to a strong bi- modality in the final black, hole mass function. 
Black hole masses below the gap can be readily understood through normal stellar processes. A 
separate and distinct population of significantly more massive black holes may result from dynamical 
processes in the core of dense globular clusters (Miller & Hamilton 2002; Giiltekin, Miller, & 
Hamilton 2004). 


4.5. Comparison with Galactic HMXBs 

Supposing that the Galaxy’s bright X-ray population is representative of low-redshift galaxies, 
we expected to find a number of sources similar to Galactic X-ray binaries in our sample. In 
our sample, we find approximately 24 sources with luminosities below our high-state ULX cutoff 
(ps. 3 x 10 39 ergs -1 ) , X-ray positions within the optical extent of their host galaxy, and no obvious 
optical counterpart. The unabsorbed luminosities for these sources range from 0.4— 2.5xl0 39 ergs -1 
(0.3 — 10 keV band) . Two of these sources were transients in the XMM data. Of the four host 
galaxies with multiple observations examined, two of these galaxies contained solely ULX sources 
in our luminosity regime (Holmberg II and NGC 5204). Each of the remaining two (NGC 253 and 
NGC 4258) had a transient source best fit by a combined blackbody and a power law. 
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This suggests an interesting diagnostic in terms of distinguishing our ULX sources from a 
normal HMXB population. In our own galaxy, most HMXBs vary on timescales of days or less and 
most of the black holes in the Milky Way are transients, though some HMXBs are indeed persis- 
tent. The figures in Kalogera et al. (2004), determined through detailed mass-transfer calculations, 
indicate that transient behavior should not be expected from a population of IMBHs. Thus, on 
average, our ULX sources should remain X-ray bright in multiple observations. Through a liter- 
ature search, we found that 37/42 of our ULX sources are well detected in ROSAT observations 
and thus are luminous for greater than 10 years and therefore are not transients. Examination of 
the long term light curves show that most of these sources vary by less than a factor of 3 over the 
timescale from ROSAT to XMM. The sources that have been above the Eddington limit in the 
Milky Way and the Magellenic clouds do so transiently, for a small fraction of the time. As best 
as we can tell, from the light curves from Einstein, ROSAT, ASCA, Chandra and XMM the ULXs 
are, rarely, transients, and are almost always ‘on’, unlike Galactic “ultra-luminous” objects. 

As a possible further diagnostic, we constructed a color-color diagram for our ULX sources. 
We adopted the colors of Done & Gierlinski (2003) in order to compare our sample with their 
sample of Galactic X-ray sources. Thus, our colors were constructed from unabsorbed model fluxes 
in four energy bands: 3-4, 4-6.4, 6. 4-9. 7, and 9.7-16 keV. The XSPEC command dummyresp was 
used to calculate a flux based on the model for the 10-16 keV range. We plot colors for a pure 
unabsorbed power law (from r = 1.5 — 3.0) and an unabsorbed MCD model ( diskbb in XSPEC mth 
kTi n = 5.0 - 0.2 eV) for comparison. Comparing our Figure 6 with Figure 8 of Done & Gierlinski 
(2003), we find that our ULX sources largely lie along the same regions as their black hole sources. 
A few ULX sources, however, lie in the region occupied by atoll and Z-sources in the plot of Done 
& Gierlinski (2003). These sources were those best fit by a Comptonization model. 


4.6. Galaxy Sample 

In this section, we examine the environment in which ULX sources reside. We investigate 
the claim that the ULX population is proportional to the host galaxy’s star formation rate (SFR) 
(Ran alii et al. 2002; Grimm et al. 2003). Towards this end, we use the far-infrared luminosity of 
the host galaxy as an indicator of the SFR. In order to compare the ULX population of a galaxy 
with the SFR we followed a similar approach to Swartz et al. (2004). We calculate the FIR flux 
from observations taken by the Infrared Astronomical Satellite. As in Swartz et al. (2004), the 
flux between 42.4 and 122.5 pm is approximated as 1.26xl0 -11 (2.58S6o+Sioo) ergcm -2 s -1 . The 
values of the flux at 60 pm (S6o) and 100 pm (Smo) were obtained from either Ho et al. (1997) or 
NED. Luminosities were calculated using the distances quoted in Table 1. We list these values in 
addition to the number of ULXs observed in individual galaxies in Table CIO. The number of ULXs 
includes both the objects we classify as high and low state ULX as well as those sources resolved 
by Chandra. 

In Figure 7, we show two plots relating the number of ULXs to L fir. It has been suggested 
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by Grimm et al. (2003) that the luminosity function in the X-ray regime from HMXBs is related to 
SFR. In our first plot, we find that the galaxies with the highest Lfir seem to have fewer ULXs than 
may be expected from the luminosity functions of Grimm et al. (2003), who present a relationship 
showing a scaling of the number of HMXBs with luminosities over a set threshold with the host 
galaxy’s SFR (see equation 7; Grimm et al. 2002). Using this relationship, we would expect that 
a galaxy with a SFR approximately equal to that of M51 4M0yr _1 from their table 1) to have 

R3 4.47 objects with luminosities greater than 10 39 ergs -1 . We find 5 objects with this luminosity in 
M51, consistent with their result. However, for NGC 4945, a galaxy with approximately the same 
L fir and therefore SFR, we find only one source with a luminosity in this range. However, we note 
that NGC 4945 is a Seyfert, implying that the L fir may primarily be caused by the AGN and 
not a direct indication of SFR. In addition to high Lfir sources with few ULX we find a number 
of sources with very small SFR but which contain a ULX. For sources with SFR < O.2M0yr -1 , 
which corresponds roughly to sources with L fir less than that of NGC 4736, we would expect 
< 0.22 sources with luminosities above 10 39 ergs _1 . However, there are a number of bright ULXs 
in galaxies with very low SFRs (for instance Holmberg II, Holmberg IX, NGC 5204, and NGC 
5408). Thus, in a direct comparison, our results do not agree with the predictions of Grimm et al. 
(2003). 

The second plot displays the average number of ULXs/galaxy, binned according to luminosity. 
This plot is extremely similar to Fig. 15 of Swartz et al. (2004) for spiral galaxies. Thus, once 
again, it seems that the connection between SFR and the ULX population in spirals is supported. 
For irregular galaxies, however, there seems to be more of a spread in the distribution. This could 
be the result of poor sampling — most of the bins contain only one galaxy. Another possibility is 
that there is no direct correlation in irregular galaxies or that the overall star formation in these 
galaxies is less ordered or clumpier. If the latter is the case, the overall SFR of the galaxy is only 
an average over a wide range of values. We shall address this issue again in the next paper in this 
series (L.M. Winter et al., submitted) where we discuss the local environments of the ULXs in our 
sample. 

In Figure 8 we plot the distribution of column densities among the ULXs. We subtracted the 
Galactic column density towards the galaxy (obtained from the nH FTOOL and listed in Table 1) 
from the values obtained through spectral fits. We note that, on average, the ULXs have large 
column densities. The typical Galactic column density along a line of sight is « 4x 10 2O cm~ 2 . 
If the ULX is located on the opposite side of its host galaxy, we might expect maximum column 
densities of « 1.2 x 10 21 cm" 2 . However, most of our sources have column densities well above this 
value. This is in agreement with the analysis of 5 ULXs by Roberts et al. (2004) and may imply, 
as they suggest, that the local environment of the ULXs contains an extra source of absorption. 
We are investigating this further, comparing the X-ray absorption column densities with HI data 
(L.M. Winter et al., submitted). 

In order to better understand the relationship between SFR and the ULX population, it is 
necessary to extend ULX studies to other wavelengths. In particular, it becomes important to 



analyze UV and IR images close to the ULX. 


5. Conclusion 

We present the results of an XMM survey of the ULX population in nearby galaxies. In this 
study, we assumed that ULXs are isotropic emitters. For our selected ULX sources (which excluded 
transient sources and supernovae) , this assumption was supported by the finding that 37/42 of our 
ULXs were found to be ‘on’ in ROSAT observations. This implies that theses sources exhibited high 
luminosities for time scales of at least 10 years, a property that is not seen in Galactic Eddington- 
limit exceeding sources (such as black hole X-ray binaries undergoing an outburst) . We also assumed 
that if some ULX sources represent a class of IMBH X-ray binaries, they would exhibit spectral 
states analogous to Galactic stellar mass black hole X-ray binaries. This is the hypothesis we set out 
to test, classifying a source as a ULX based on (1) spectral form, (2) luminosity, and (3) coincidence 
of the X-ray source within the optical host galaxy. Due to the quality of spectra available for these 
distant X-ray sources, our classification of spectral form is really a first approximation describing 
the basic curvature of the spectrum. 

Through this study, we have found that there exists a population of objects whose X-ray 
spectral properties closely match the low/hard state spectra of Galactic black holes, but whose 
luminosities lie in the range of Lj, 0 ; w 2 X 10 38 — lx 10 4 ° ergs -1 . In the Milky Way, black holes 
with these spectral properties radiate at only w 0.05 of the Eddington limit. If this is also true for 
this population, it indirectly implies that these objects have a mass greater than « 30 Mg ranging 
up to 1500 Mq and thus should be IMBHs. The existence of such objects was “predicted” on 
the basis that the ULXs previously studied shared the X-ray spectral characteristics of high-state 
Galactic black holes; namely, an X-ray spectrum best fit by a combined blackbody and a power 
law (Miller et al. 2003), but with much higher luminosities. If these objects are high-state IMBHs, 
the corresponding low-state objects should also exist. 

Our survey has also uncovered a large population of objects whose X-ray spectra are well 
modeled by the canonical description of Galactic black holes in the high-state (thermal dominated), 
a black hole with a steep power law, but whose bolometric luminosities exceed 2 x 10' 39 ergs -1 , 
ranging up to 10 41 ' 5 ergs -1 and whose blackbody temperatures are less than 0.3 keV. If these 
objects are radiating at « 1/2 the Eddington limit like their Milky Way counterparts their implied 
masses are from 30 — 3000 M©, a range very similar to that implied by the low-state objects. 
Using the M” 1 / 4 scaling of mass to temperature, the observed spectral temperatures give masses 
of 500 — 10 4 M© a considerably larger value. In general agreement with the expectations of the 
IMBH hypothesis, the objects with high-state spectra are more luminous than those with low-state 
spectra. We note that these results have required the high signal to noise of XMM in order to 
discern the spectrum of these objects. Many of these objects have also been observed by Chandra 
and their spectra have been well-fitted by simple power laws. 
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In addition to classification of the sources, we investigated some of the properties of the ULX 
sources. We found a gap in the temperature distribution of high/soft state ULXs. This gap may 
indicate a gap in mass distribution, which may provide clues to the nature of ULXs. We also 
found that our ULXs are persistent sources (not transients) which occupy regions on the color- 
color diagram of Done & Gierlinski (2003) also occupied by Galactic black hole sources. Lastly, the 
existence of a substantial population of ULXs in nearby dwarf and other low star formation rate 
galaxies argues that (in agreement with Ptak & Colbert (2004); Swartz et al. (2004)) there is more 
than one source term for the origin of ULXs, with at least some of them not being associated with 
recent star formation, at least statistically. 

We conclude, from an X-ray spectral and luminosity point of view, that our data are consistent 
with many of these objects having the properties expected of an IMBH population. However, we 
also find two other populations of objects, those whose blackbody temperature and luminosity 
correspond to that of stellar mass black holes with kT s-s 1 keV and log Lx less than 2 x 10 39 erg s -1 
and a small population of objects whose X-ray spectra and luminosities are consistent with that 
of stellar mass black holes in the very high state. Thus, ULX selected purely on the basis of 
0.3 - lOkeV X-ray luminosities are a composite class with 1/4 being “normal” stellar mass black 
holes and the rest being consistent with a population of IMBHs. 

In a follow-up paper we will discuss the environments of these objects as revealed by XMM 
OM UV imaging and the implications this has for the origin of ULXs. 


L.W. gratefully acknowledges Kip Kuntz and M. Coleman Miller for helpful discussions. We 
would also like to acknowledge R. Narayan for asking the question “are there any low-state ULXs?” 
at the Kyoto Black Hole meeting. 


REFERENCES 

Anderson, S. et al. 2003, AJ, 126, 2209 

Bauer, F.E., Brandt, W.N., & Lehmer, B. 2003, AJ, 126, 2797 

Begelman, M. C. 2002, ApJ, 568, L97 

Belloni, T., et al. 1997, ApJ, 488, L 109 

Brandt, W.N., Iwasawa, K., & Reynolds, C.S. 1996, MNRAS, 281, L41 
Chen W., Schrader C. R., Livid M., 1997, ApJ, 491, 312 
Colbert, E.J.M. & Ptak, A.F. 2002, ApJS, 143, 25 
Done, C. & Gierlinski, M. 2003, MNRAS, 342, 1041 



- 21 - 


Ebisawa, Iv., et al. 2003, ApJ, 597, 780 

Eracleous, M., Shields, J.C., Chartas, G., & Moran, E.C. 2002, ApJ, 565, 108 
Freedman, W. L., et al. 1994, ApJ, 427, 628 

Freeman, K. C., Karlsson, B., Lynga, G., Burrell, J. F., van Woerden, H., & Goss, W. M. 1977, 
A&A, 55, 445 

Gierliriski, M. & Done, C. 2004, MNRAS, 349, L7 

Grimm, H.-J., Gilfanov, M., & Sunyaev, R. 2003, MNRAS, 339, 793 

Guainazzi, M., et al. 2000, A&A, 356, 463 

Giiltekin K., Miller M.C., Hamilton D.P., 2004, ApJ, 616, 221 

Gutierrez, C.M. & Lopez-Corredoira, M. 2005, preprint (astroph/0502290) 

Hartwell, J., et al. 2004, MNRAS, 348, 406 
Heger A., Woosley S.E., 2002, ApJ, 567, 532 

Ho, L.C., Filippenko, A.V., & Sargent, W.L.W. 1997, ApJS, 112, 315 

Humphrey, P.J., Fabbiano, G., Elvis, M., Church, M.J., & Balucinska-Church, M. 2003, MNRAS, 
' 344, 134 

Jenkins, L.P., Roberts, T.P., Warwick, R.S., Kilgard, R.E., & Ward, M.J. 2004, MNRAS, 349, 404 
Jones, C. 1977, ApJ, 214, 856 

Jonker, P.G. & Nelemans, G. 2004, MNRAS, 354, 355 

Kaaret, P., Corebel, S., Prestwich, A.H., & Zezas, A. 2003, Science, 299, 365 

Kalogera, V., Henninger, M., Ivanova, N., & King, A.R. 2004, ApJ, 603, L41 

Karachentsev, I.D., Sharina, M.E., Dolphin, A.E., et al. 2002, A&A, 385, 21 

Kelson, D. D. 1996, ApJ, 463, 26 

King, A.R. et al. 2001, ApJ, 552, L109 

Kong, A.K.H. 2003, MNRAS, 346, 265 

Kong, A.K.H. & Di Stefano, R. 2003, ApJ, 590, 13 

Kording, E., Falcke, H., & Markoff, S- 2002, A&A, 382, L13 



22 - 


Kubota, A. & Makishima, K. 2004, ApJ, 601, 428 
Kubota, A. et al. 2001, ApJ, 560, L147 
Kuntz, K. et al. 2005, ApJ, 620, L31 
Lehmann, I. et al. 2005, A&A, 431, 847 

Lira, P., Lawrence, A., & Johnson, 11. A. 2000, MNRAS, 319, 17 
Liu, J. F., Bregman, J. N., & Seitzer, P. 2004, ApJ, 602, 249 
Madau P., Rees M.J., 2001, ApJ, 551, L27 
Magdziarz, P. & Zdziarski, A. A. 1995, MNRAS, 273, 837 
Makishima, K., et al. 2000, ApJ, 535, 632 

Martin, C.L., Kobulnicky, H.A., fc Heckman, T.M. 2002, ApJ, 574, 663 
McClintock, J.E. & Remillard, R.A. 2004, preprint (astroph/0-306213) 

Miller, J. M., Fabbiano, G., Miller, M. C., & Fabian, A. C. 2003, ApJ, 585, L37 

Miller, J.M., Fabian, A.C., & Miller, M.C. 2004, ApJ, 614, L117 

Miller, M. C. & Colbert, E. J. M. 2004, IJMPD, 13, 1 

Miller M.C. fe Hamilton D.P., 2002, MNRAS, 330, 232 

Miyamoto, S. et al. 1991, ApJ, 383, 784 

Pietsch, W. et al. 2001, A&A, 365, 174 

Pietsch, W., Haberl, F., & Vogler, A. 2003, A&A, 402, 457 

Ptak, A. & Colbert, E. 2004, ApJ, 606, 291 

Pakull, M.W. & Mirioni, L. 2003, RMXAA, 15, 197 

Porquet, D., Reeves, J.N., O’Brien, P., & Brinkmann, W. 2004, A&A, 422, 85 
Radecke, H.D. 1997, A&A, 319, 18 
Reynolds, C.S., et al. 1997, MNRAS, 286, 349 
Roberts, T.P. & Warwick, R.S. 2000, MNRAS, 315, 98 

Roberts, T.P., Warwick, R.S., Ward, M.J., & Goad, M.R. 2004, MNRAS, 349, 1193 
Roberts, T.P. et al. 2005, MNRAS, 357, 1363 



Ross, R.R. & Fabian, A.C. 2005, MNRAS, 358, 211 

Sakai, S., Madore, B.F., & Freedman, W.L. 1996, ApJ, 461, 713 

Sakai, S. & Madore, B.F. 1999, ApJ, 526,599 

Schlegel, E.M., Petre, R., Colbert, & Miller, S. 2000, AJ, 120, 2373 

Schlegel, E.M. et al. 2004, ApJ, 603, 644 

Shapley, A., Fabbiano, G., & Eskridge, P.B. 2001, ApJS, 137,139 
Schlegel, E.M., Petre, R., Colbert, E.J.M., & Miller, S. 2000, AJ, 120, 2373 
Smith, D.A. & Wilson, A.S. 2001, ApJ, 557, 180 
Stocke, J.T. et al. 1983, ApJ, 273, 458 

Summers, L.K., Stevens, I.R., Strickland, D.K., & Heckman, T.M. 2003, MNRAS, 342, 690 

Summers, L.K., Stevens, I.R., Strickland, D.K., & Heckman, T.M. 2004, MNRAS, 351, 1 

Sunyaev, R.A. & Titarchuk, L.G. 1980, A&A, 86, 121 

Swartz, D.A., Ghosh, K.K., Tennant, A.F., & Wu, K. 2004, ApJS, 154, 519 

Terashima, Y. & Wilson, A.S. 2004, ApJ, 601, 735 

Tolstoy, E., Saha, A., Hoessel, J.G., & McQuade, K. 1995, AJ, 110, 1640 

Tosi, M., Bellazzini, M., Aloisi, A., Greggio, L., Leitherer, C., & Montegriffo, P. 2001, AJ, 122, 
1271 

Tully, R.B. 1988, Nearby Galaxies Catalog. Cambridge Univ. Press, Cambridge 
Voges, W. et al. 1999, A&A, 349, 389 
Vogler, A. & Pietsch, W. 1996, A&A, 311, 35 
Wang, D. et al. 2004, ApJ, 609, 113 

Zdziarski, A. A .& Gierlinski, M. 2004, Progress of Theoretical Physics Supplement, 155, 99 
Zimmermann, H.U. & Aschenbach, B. 2003, A&A, 406, 969 


This preprint was prepared with the AAS I£TgX macros v5.2. 



- 24 - 


Table 1. XMM-Newton Galaxy Observations 


Galaxy 

Type a 


distance 0 

ref 

obs id d 

duration (s) 

comments 

NGC247 

SAB(s)d 

1.54 

3.09 


0110990301 

14536 

- 

NGC253 

SAB(s)c;HII 

1.40 

3.73 


0110900101,0152020101 

30711,110591 

Starburst 

NGC300 

SA(s)d 

3.11 

2.56 


0112800101 

43967 

- 

NGC625 

SB(s)m? sp; HII 

2.15 

2.62 


0085100101 

26288 

- 

NGC1313 

SB(s)d; HII 

4.0 

4.17 


0106860101 

41310 

- 

IC0342 

SAB(rs)cd; HII 

30.3 

3.9 

1 

0093640901 

11217 

- 

NGC1569 

IBm 

21.7 

1.6 

1 

0112290801 

15582 

Starburst 

NGC1705 

SA0- pec; HII 

3.9 

5.1 

2 

0148650101 

58926 

Starburst 

MRK 71 

BCD; HII 

3.9 

3.4 

3 

0141150201 

45919 

galaxy pair 

NGC2403 

SAB(s)cd; HII 

4.15 

3.56 


0150651201 

11415 

- 

Holmberg II 

Im 

3.42 

2.70 


0112520701,0112520901 

13528, 6860 

- 

Holmberg I 

IAB(s)m 

3.49 

3.6 

4 

0026340101 

26280 

- 

M81 

SA(s)ab;LINER 

4.12 

3.6 

4 

0111800101 

127913 

Hoi IX also in field of view 

M82 

10; HII 

4.14 

3.9 

5 

0112290201 

29387 

Starburst 

Holmberg IX 

Im 

4.0 

3.6 

4 

0112521001 

10350 

MSI also in field of view 

Sextans A 

IBm 

3.85 

1.4 

6 

0026340201 

21618 

- 

IC 2574 

SAB(s)m 

2.29 

3.6 

7 

0026340301 

24263 

bursting star-formation 

NGC 4214 

IAB(s)m; HII 

1.49 

2.7 


0035940201 

14744 

- 

NGC 4258 

SAB(s)bc;LINER 

1.2 

7.2 


0059140901, 0110920101 

16146, 21895 

- 

NGC4395 

SA(s)m;LINER 

1.33 

4 


0112521901 

15842 

- 

NGC4449 

IBm; HII 

1.39 

3.08 


0112521701 

15522 

- 

NGC4490 

SB(s)d 

1.78 

7.8 

1 

0112280201 

17754 

interacting with NGC4485 

NGC4631 

SB(s)d 

1.28 

7.5 

1 

0110900201 

53850 

- 

NGC4736 

(R) S A (r ) ab ;LINER 

1.43 

4.3 

1 

0094360601 

23461 

- 

NGC4945 

SB(s)cd; Sy2 

15.9 

3.1 


0112310301 

23062 

- 

NGC 5204 

SA(s)m; HII 

1.42 

4.8 

1 

0142770101, 0142770301 

19205, 16387 

- 

M51 

Sc; Sy2 

1.55 

7.2 

1 

.0112840201 

20924 

Galaxy pair 

M83 

SAB(s)c;HII 

3.94 

6.2 


0110910201 

30627 

Starburst 

NGC5253 

Im pec;HII 

3.77 

3.2 

1 

0035940301 

47216 

Starburst 

M101 

SAB(rs)cd 

1.17 

7.4 

8 

0104260101 

43019 

- 

NGC5408 

IB(s)m; HII 

5.73 

4.8 

9 

0112290601 

7757 

- 

Circinus 

SA(s)b; Sy2 

57.8 

4 

10 

0111240101 

110496 

- 


a from the NAS A/IP AC Extragalactic Database (NED) 

b column density in units of 10 20 cm" 2 , obtained from the web version of the nH FTOOL 
c distance in Mpc (if no reference is given, obtained from the distance modulus given in LEDA) 
d XMM-Newton observation ids for the data examined in this survey 


References. — (l)Tully 1988; (2) Tosi et al. 2001; (3) Tolstoy et al. 1995; (4) Freedman et al. 1994; (5) Sakai Si Madore 1999; (6) Sakai, 
Madore, & Freedman 1996; (7)Shapley, Fabbiano, &: Eskridge 2001; (8) Kelson 1996; (9) Karachentsev et ah 2002; (10) Freeman et al. 1977. 
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Table 2. XMM-Newton best fit: single component spectral fits 


Powerlaw Bremsstrahlung 


Source 


r 

X 2 /dof 


kT (keV) 

X 2 /dof 



NGC247 XMM2 

i-4±i;t 

9 oq+1-02 

^y_0.57 

47 . 7/54 

< 0.65 

9 c:c+6-90 

^• 00 — 1.61 

9 1 9+0-52 
4i*l^_o.37 

48.8/54 

0.33 

0.38 

NGC 253 XMM2 (obs 1) 


o si+°- 18 
■ *^> A — 0.17 

69.1/74 

o-5l8;! 

74.7/74 

0.52 

0.87 

NGC300 XMM4 d 

2.5 

9.07 

90.6/45 

0-27 

0.14 

117.6/45 

- 

- 

NGC 1313 XMM4 

1.86±°-° 

i-s+255 

141.7/149 

1 2+°' 3 
A * Z -0.3 

6.62+! ;| 8 

140.1/149 

0.33 

0.69 

IC0342 XMMl 

5.8i° 0 ;i 

1J68.+8;2S 

159.5/185 

4-9+85 

io-5l!;| 

160/185 

3.5 

6.37 

IC.0342 XMM2 

23.9l|;“ 

1 85+ 0 * 22 
4-©i>_ 0 .20 

77.5/85 

21+15 

8.5±|g 

74.9/85 

4.64 

8.44 

IC0342 XMM4 

5-3li;I 

2-02i°j° 

64/58 

4 o+O-OO 
4,Z -0.85 

4-441!;°° 

56.9/58 

0.69 

1.26 

MRK71 XMM1 

0 47+ 0 - 30 
-0.32 

i-69l°;“ 

55.3/54 

0 04^°* 26 
■ — 0.04 

7 -98l!;!8 

59/56 

0.19 

0.27 

NGC2403 XMM4 

1 y+O.S 
A *'-0.7 

i-89l°;|° 

62.3/71 

i-i±g5 

4.591!;! 

62.3/71 

0.31 

0.48 

HolmlIXMMI (obs 2) 

1 5+°- 2 
A '°-0.2 

3.09l°;“ 

266.7/252 

°-3il8;i! 

1 13+°- li 
A * AC5 -o.ii 

309.4/252 

3.5 

3.1 

Holm I XMM2 

0.35 e 

2-i3l°;l° 

39.2/45 

0.35 e 

2 0 + 0- 57 

■ * y — 0 .48 

68.8/45 

0.10 

0.16 

Holm I XMM3 

0.35 e 

9 QC-+ 0.1:9 

34.4/32 

0.35 s 

9 no+0.85 

42.1/32 

0.12 

0.19 

IC2574 XMMl 

i O+0.40 

1 Q 7+ 0-0.7 
A ' y ' -0,10 

120.9/103 

■0.69±g;g 

4 4+0.89 
^* A — 0.67 

107.5/103 

0.35 

0.47 

IC2574 XMM2 

0-4tH 

9 9 + 0.21 
^-0.09 

45.7/51 

0.229 e 

7 71 +0-33 
* * ■ A — 0.27 

65.4/52 

0.22 

0.34 

XC2574 XMM3 

O-IsIo'm ’ 

9 4Q+0* 27 

40.3/49 

0.229 e 

0.97i°!° 

76.5/49 

0.22 

0.34 

NGC4214 XMMl 

i r +0.52 
A * A — 0.47 

1 07 +O .26 

i ’ tSY - 0 . 2 l 

41.9/38 

0-541°;!! 

4.86t!;f! 

44.5/38 

0.25 

0.22 

NGC4258 XMM2 (obs 2) 

6 7 +2 ' 6 
°' f -1.5 

2 49+° - 44 

83.6/57 

4.8l°;! 

2.6il!; 2 7 | 

85.5/57 

0.30 

1.9 

NGC4258 XMM3 

1 4+0-69 
1 ^_0.64 

9 OO + 0.34 
— 0.24 

38.9/37 

0 49+ 0 * 44 
u *"*y_o.38 

2.481!;°! 

41.3/37 

0.20 

1.2 


3.8 

1.82 

4/11 

2.7 

7.14 

5/11 

0.077 

0.48 

NGC4258 XMM4 

o.esl®;^ 

1 07 + 0-22 
A * y ( -0.19 

41.1/48 

u-Ub_ 0 05 

4.071!;° 

45.2/48 

0.39 

2.4 


1 n+0-78 
A ‘ J -0.60 

9 94+0- 2 0 
^•^4-0.24 

77.03/77 

q 9+0.6 
U ’ y -0.4 

2.821!;! 

77.8/77 

0.33 

2.0 

NGC4395 XMM2 

0-33l°;| 

9 7C+0.45 
z - ' ° —0.33 

38.6/36 

0.133® 

0.79+811 

52/37 

0.15 

0.28 

NGC4395 XMM4 

0.3l®j 

2-08l°;|° 

16/25 

0.133 s 

1 Q7 + 0-99 

A * y < —0.60 

21.6/26 

0.15 

0.28 

NGC4449 XMM2 

l-5l°o:! 

9 81 +0 16 

103.5/112 

o.25t°;! 

i fijr + 0.22 
A * ■ — 0.21 

112.1/112 

0.29 

0.33 

NGC4490 XMM4 

10,2+ff 

2.09l°;^ 

51.6/50 

8.3li;l 

4-75+!;!! 

50.3/50 

0.84 

6.1 

NGC4490 XMM5 

o q+0.94 
^■^-0.81 

9 91+0-22 
^ ■ A — 0.20 

60.1/65 

2-5lg;J| 

3.08l°;|! 

61.6/65 

0.41 

2.98 

NGC4631 XMM4 d 

7.8 

9.50 

261.5/74 

2.9 

0,17 

207.8/74 

- 


NGC4631 XMM5 f 

1.3 

1.03 

641.8/153 

1.3 

199 

659/153 


- 

NGC4945 XMM3 

3-3±|;| 

1 8940-12 
A 0.20 

30.1/30 

9 c+0-63 
Z ‘°-0.90 

6.07+ 4 '* 9 

30.3/30 

0.38 

0.43 

NGC5204 XMM2 

0-89+°j|, 

1 QR+0-25 

x.ye_ 0 _2o 

42.37/42 

0.23±g;| 2 

4 0 5 +l ; M 

42.2/42 

0.15 

0.41 


0.7B±S'.45 

l-63lo'.i7 

41.4/47 

0.42l°;|° 

7.821!;!! 

39.4/47 

0.25 

0.69 

M51 XMM3 

°-6-o'.40 

1-S6±g;53 

63.2/72 

■0.05+8;g a 

5.22l 2 -j° 

69.2/72 

0.18 

1.1 

M51 XMM4 

■0-4+SS 

1 cr+0.08 

l * o ^_0.13 

34.8/37 

0.0118;!! 

111 +0-32 
AA - A -0-25 

34.8/37 

0,16 

0.99 


a total column density in units of 10 21 cm -2 

b unabsorbed flux in the 0.3-10 keV band in units of 10 -12 ergcm“ 2 s - * 1 

c unabsorbed luminosity in the 0.3-10 keV band, using the distances quoted in Table 1, in units of 10 39 ergs' 
d see appendix; super-soft X-ray source best fit by single-component biackbody 
e absorb tion column density fixed to the galactic column density found in Table 1 
f source is best fit by a combined power law and vapec model; see appendix 
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Table 3. XMM-Newton best fit: two-component blackbody and power law spectral fits 


Source 


kT (keV) 

r 

X 2 / dof 

Ax 2b 

Fx c 


NGC247 XMMl 

4 1 + 1 ’ 9 

0.l2±lf 2 

4.181!;!! 

86.5/93 

25.7 

6.2 

7.1 

NGC253XMM1 

2 <7+0 A 

o-BPiSiS 

i-wjoj I 

225.9/230 

36.7 

2.7 

4.5 


7 3 + 1,1 

1 14+0-07 
A * 4 ^- 0.10 

2 54+0 * 27 

567/580 

44.6 

3.4 

5.7 

NGC253 XMM 2 (obs 2 ) 

2 0+°- 3 
■ * ■ — 0.2 

0 71+0.1° 
U * / 4 - 0.10 

O 14+0.05 
Z * A ^-0.08 

460.3/498 

47.1 

1.6 

2.7 

NGC253 XMM3 

3 1 + 4,8 
- 4 — 0.5 

0 75+0-13 
U * ' - 0.10 

2 47 + 2 - 99 

68.5/82 

23.4 

0.60 

1.0 


3- 2 ±o:s 

0 67 4 * 0 ' 13 
u * b '-0.09 

2 07+ 0 - 14 
*- Uf -0.20 

347.4/407 

34.4 

0.80 

1.3 

NGC253 XMM4 

20 i ^ 8 

0 - 11 Io’.O3 

2 51 +0.49 
^* bA — 0.30 

66.7/57 

6.9 

15 

25 


4 5 "^ 1,2 

0.09l°;^ 

9 on+0.27 

309.3/291 

12.1 

1.4 

2.2 

NGC253 XMM5 

i c + 7.2 
4 ' b — 1 -5 

0.96i°;l| 

2.431!;°® 

26.5/23 

5.3 

0.26 

0.43 



°- 16 -0'.03 

1 Q5+0.-1 4 
1,95 x 0 11 

223.7/296 

60.1 

1.4 

2.2 

NGC253 XMM 6 

6-3±?;i 


9 26 4 * 0 ' 48 

tRJ? 

417.9/407 

17.1 

1.9 

3-1 

NGC253 XMM7 

6.3±!;! 

0.69l°;“ 

9 40 + 0 - 4 1 

335.8/339 

21.2 

1.1 

1.8 

NGC300 XMMl 

i 7 + 0.20 

0.98+S:lo 

3 41+ 0 * 06 

443.7/420 

26.1 

1.3 

1.0 

NGC300 XMM 2 

3-8ll ; 7 4 

0.091°;°! 

2.87l° 0 ; 3 3 4 8 

102.6/97 

31.34 

1.1 

0.86 

NGC300 XMM3 


°-04i85J 

1.981!;! 

87.7/79 

14,2 

1.2 

0.93 

NGC300 XMM 6 

9 o+ 2.6 

^-1.3 

0.84±°;^ 

4-9±o !? 7 

34.6/35 

13 

0.27 

0.20 

NG Cl 313 XMMl 

3 0+*- 2 
8 - U -0.9 

0 13+ 0 - 03 
U * A ' 3 - 0.02 

1.751°;!! 

194.1/201 

35.4 

0.64 

1.3 

NGC1313 XMM 2 

o i+ 0 - 4 

0 1 6+0-04 
U * lb - 0.02 

9 97 + 0.10 

TH 2 

425.2/419 

38.9 

2.0 

4.2 

NGC1313 XMM3 

0 < 7+0 .8 
b - 4 -o .6 

o-ulSSl 

2.761!;!! 

441.7/424 

336.6 

10 

22 

IC0342 XMM3 

9 7 + 1 - 8 

0.09l°;°! 

2.691!;!® 

129.5/107 

56.3 

31 

56.4 

NGC1705 XMMl 

A 9Q+0-39 

1.01±o'.29 

2. 31 ±|-|1 

53/85 

8.9 

0.10 

0.44 

NGC1705 XMM3 

<1.44 

i,orl°;!° 

9 2 Q+ 0.70 

69.8/65 

11.1 

0.15 

0.48 

NGC2403 XMMl 

2.3±i:f 

o. 66 l°;!° 

9 1 q+0.41 
^* 4 ° — 0.59 

81.4/79 

10.8 

1.99 

3.1 

NGC2403 XMM 2 

1 0 + 0.8 
4 ‘ 8 - 0.6 

0.62l°;!! 

1 95 4-0,26 
l- y i>_0.42 

163.1/151 

16.4 

1.0 

1.6 

NGC2403 XMM3 

1 7 + 1 * 1 
4 *'- 0.8 

0 74+0-23 

9 1 5 +O .66 

^■la_0.40 

84.2/105 

8.4 

0.64 

1.1 

Holm II XMMl (obs 1 ) 

1 6 +0 1 

0.i4l°;°! 

9 35 + 0.05 

997.5/976 

136.7 

12 

10 

M81 XMMl 

o O+0.17 

^-o.os 

0.90l°;°| 

9 52+0-03 
4 * D4 -0.04 

1273.7/1243 

533.1 

4.5 

7.0 


3 5+ 0 ' 4 
8 * b -0.6 

1 i3+ 0 -i 3 
A * A ^— 0.14 

9 34+0- 29 

203.5/204 

21.4 

4.8 

7.4 

MSI XMM 2 

7 4 + 0 * 5 

a -t +0-004 
U - A -0.004 

9 07 + 0.16 
tHX 

833.9/616 

524.3 

13 

22 

MSI XMM4 

1 l 4 * 1 * 6 
1 , 1 —IT) 

9 51 +I-II 
j 2 * b 4 -0.73 

9 31 + 1 - 22 

48.9/50 

28.2 * 

0.43 

0.70 

MSI XMM5 

n 1t r+0-69 
u * 4 b -o.i3 

0 62 4 * 0 19 
U * b 4 -D.ll 

1 26+ 0 * 22 

89/80 

8.5 

0.38 

0.59 

Holm IX XMMl 

9 I + 0.2 
- *' l -0.2 

0 ^ 7+0-02 
U * A ' — 0.02 

1 79+0-04 
4 - 4 — 0.03 

866.6/878 

134.3 

10 

16 

NGC4258 XMMl 

0.38±!;!® 

o.54l °;! 7 


91.1/76 

10.3 

0.34 

2.1 

NGC4258 XMM 2 (obs 1 ) 

1 9 + 2 * 4 
A -* y -0 A 

0 - 78 to '.!3 

2 . 02+!'! 5 

73.4/61 

24.1 

0.31 

1.9 

NGC4395 XMMl 

9 0 + 0.08 
Z * U -0.07 

O.i4l°o ;° 0 2 2 

3 44 + 0 -' 54 

168.2/154 

26.9 

1.4 

2.7 

NGC4449 XMMl 

8 - 7 x?i 

0 ' 19 zL°! 

9 91 +0-33 

rS-if 

111.2/116 

4.3 

1.56 

1.8 

NGC4449 XMM3 

3-3i:l 

n 1 5+0.03 
U * 15 -0.03 

2.52l!-j® 

119.9/87 

34.1 

1.1 

1.3 

NGC4490 XMMl 

r 0 + 2.96 
3 *s_2.96 

a 77+0.08 
^ ' —0.095 

2 - 89 ig|I 

66.5/63 

35 

0.88 

6.4 

NGC4631 XMMl 

O+0.9 

^-O.S 

o.i 2 i°;°! 

9 1 9+O.O3 
4 - 4/ t 0.02 

371.3/345 

12.1 

0.96 

6.5 

NGC4631 XMM 2 

2 3 + 4 * 4 

0 . 181 °;°® 

1-801!;!! 

107.4/97 

12.1 

0.25 

1.7 

NGC4631 XMM3 

1 l 4 ' 11 
L 1 - 0.8 

l.oilol 2 

2.451J 62 

127.1/96 

18.9 

0.15 

1.0 

NGC4945 XMMl 

3 5 + 2-4 

0 77+0.27 
■'* '- 0.10 

1 60 4-0 * 40 
4 0.31 

96.1/120 

20 

0.59 

0.68 

NGC4945 XMM 2 

3 2 4 ’ 1 ' 1 
0.7 

1 i 5+0*28 

4 * 4 b 70 .33 

1 cn+0-20 

4 - 80 -0.30 

105.8/113 

8.7 

0.66 

0.76 

NGC4945 XMM4 


0 61 +0-10 
U * b4 ~0.10 

9 Q 2 + 1-00 

58.4/60 

17.1 

0.38 

0.44 

NGC5204 XMMl 

0 - 66 ±°|| 

n i 0 + 0.02 
U - 4b — 0.03 

1 Q 9 + 0-12 

4 - ■ 4 — 0.06 

543.0/559 

49.1 

1.98 

5.6 


7 1 +0.08 
A ‘ 4 -0.14 

0 16+ 0,02 
b,4b — 0.02 

2 02+0-12 

461.4/496 

71.6 

2.92 

8.0 



Table 3 — Continued 


Source 

n 

kT (keV) 

r 

X 2 /dof 

Ax 2b 

Fx c 

Lx d 

M51 XMMl 

0.95lJ;i“ 

O.lslS'.os 

9 ir+ 0.42 

97/80 

13.4 

0.31 

1.9 

M51 XMM5 

10-4^;^ 

0.078^00, 

2 26"^°' 26 
2.zo_ 0 25 

59.8/70 

196.2 

220 

1900 

M83 XMM4 

i.77tf : ? 7 

09 

1 «i+ 0 - 96 

84.8/89 

6.6 

0.2 

0.92 

M101 XMMl 

0 92+ 0 ' 12 

n 91 + 0.03 
U.^l_ 0.04 

1-42 ±lit 

249.9/231 

53.1 

0.45 

2.9 

M101 XMM2 

1 6+°- 46 
A ‘ b -0.21 

0 76+ 0 - 14 

i.ss±g;?* 

251.6/261 

37.2 

0.7 

4.6 

NGC5408 XMMl 

o-9±S:i6 

0 14+ 001 

9 71 + 0-16 
z ' /J --0.20 

316.4/337 

80.4 

3.97 

10.9 

CIRCINUS XMMl 


0 10+ 001 

9 30+° - 08 
^.OU_o 08 

749.4/861 

13.5 

12 

23 

CIRCINUS XMM2 

11.2 

n co+ 0.03 

u.oo_o.03 

4-7i+“;l| 

438.5/430 

79.4 

5.6 

10.7 

CIRCINUS XMM3 

13.5+ti 

O-h7_ 0 Q8 

6.77%? 

269.3/260 

15.9 

7.6 

14.5 


Hotal column density in units of 10 21 cm -2 

b improvement in x 2 over the single-component power law model 
c unabsorbed flux in the 0.3-10 keV band in units of 10“ 12 erg cm" 2 s~* 

d unabsorbed luminosity in the 0.3-10 keV band, using the distances quoted in Table 1, in units 
of 10 39 erg s— 1 
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Table 4. XMM- Newton two-component blackbody and power law spectral fits for sources with 

large uncertainty 


Source 

n^ a 

kT (keV) 

r 

X 2 / dof 

Ax 2b 

Fx c 

Lx d 

NGC300 XMM5 

0-41 ISIS 

■I■■UO_ 0 20 

2.78+gS 

46.6/53 

7.6 

0.17 

0.13 

NGC1705 XMM2 


0 23+° 10 

1.60+gS? 

85.5/74 

6.5 

0.09 

0.27 

Holm I XMMl 

0 4 +0 - 5 
U * 4 -0.3 

1 Q 7 + 0-66 

A - y * w 0.89 

2 4g + 0 * 44 

Z * 4t5 -0.40 

97.4/93 

5.4 

0.6 

0.93 

M81 XMM3 

3 7+ 2 - 4 
■ ' —2.1 

nil +0-05 

0.1 J._o.02 

1 fiQ+0- 27 

1*UU_0 3 3 

77.1/78 

4.25 

0.53 

0.82 

Sextans A XMMl 

0 4 +°‘ 7 
U ' 4 - 0.1 

i ne+2.3 

l.UD_ 0 Q7 

2 6+ 0 ' 3 
^• b ~0.2 

269.1/271 

2.3 

0.60 

0.14 

NGC4214 XMM 2 

1 R+l - 9 
X * 8 - 0.6 

D 31 +0-56 
U - 1 — 0.21 

0 qc+l-Si 
1.05 

46.4/44 

4.5 

0.4 

0.35 

NGC4395 XMM3 

o s+°- 9 

u - & — 0.3 

1.10+555 

2.66iJ;?| 

52/56 

3.9 

0.29 

0.56 

NGC4490 XMM 2 

4 4+ 1 ' 9 

0.60 +S;S 

9 1 O+0.50 
z * lo —.0.70 

42.4/54 

7.1 

0.65 

4.7 

NGC4490 XMM3 

13+55 

o.09i“;®i 

3 91 +0-52 
c5 * zl -0.17 

72.1/78 

4.6 

12 

87.4 

NGC4736 XMMl 

6.3l|;? 

0.08t°;°l 

2.4i+5:l^ 

54.9/51 

7.9 

8.1 

17.9 

M51 XMM2 

1-3±o.5 

o.26+®;55 

l.80±S;|| 

70.7/68 

4.5 

0.36 

3.0 

M51 XMM6 

8-2.115 

0.O8+S5S 

0 n +0.37 
cS,u -0.43 

36.9/41 

4.07 

5.6 

35 

M51 XMM7 

2;8+*1 

n in +°- 03 

U.1U_ 0 Q3 

1 Q7+0-43 
A *v * —0.30 

31.7/29 

6.1 

0.26 

1.6 

MS3 XMMl 

1 g + 0.48 
1*0-0.45 

o .7415:1! 

2 SR "^°' 60 
z - O6 - 0 .24 

177.7/209 

4.7 

0.63 

2.5 

M101 XMM3 


0.6Stg;|5 

9 qq+015 

z.yo_ 0 26 

145.5/131 

3.4 

0.56 

3.7 

M101 XMM4 

1 R+0.17 

o C 4 +O.H 
V.t>4 _ 0 Q7 

9 OO+0-12 

z - zz -0.08 

158.2/138 

7.5 

0.34 

2.2 

Ml 01 XMM5 

1 3 +1 - 2 

n 10+0*05 

O.l»_ 0 .o 6 

1 ntr+0.3 

+ -J&-0.22 

45.1/44 

2.8 

0.13 

0.85 


a total column density in units of 10 21 cm” 2 

improvement in x 2 over the single-component power law model 
c unabsorbed flux in the 0.3-10 keV band in units of 10” 12 erg cm” 2 s” 1 

d unabsorbed luminosity in the 0.3-10 keV band, using the distances quoted in Table 1, in units 
of 10 39 erg s” 1 
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Table 5. Bolometric Luminosities of ULX sources 


Source 

^cutof fpl* 

Uol b 

M Bdd c 

class d 

NGC247 XMMl 

13.4258 

7.07734 

54 

HS ULX 

NGC253 XMMl 

9.31469 

2.44574 

19 

HS ULX 

NGC253 XMM2 

4.3701 

2.15292 

17 

HS ULX 

NGC253 XMM6 

5.05828 

3.92514 

30 

HS ULX 

NGC1313 XMM3 

37.0364 

27.9692 

215 

HS ULX 

NGC1313 XMM4 

1.50345 


116 

LS IMBH cand. 

IC0342 XMMl 

14.1215 


1086 

LS ULX 

IC0342 XMM2 

19.8129 


1524 

LS ULX 

IC0342 XMM3 

114.015 

95.4068 

734 

HS ULX 

MRK71 XMMl 

0.2993 


23 

LS IMBH cand. 

NGC2403 XMMl 

4.1497 

2.14873 

17 

HS ULX 

NGC2403 XMM4 

0.57068 


44 

LS IMBH cand. 

Holmberg II XMMl 

0.88906 


68 

LS IMBH cand. 


16.8335 

11.4543 

88 

HS ULX 

Holmberg I XMM2 

10.5158 


809 

LS ULX 

MSI XMMl 

15.7004 

3.17932 

24 

HS ULX 

Holmberg IX XMMl 

31.0582 

28.1445 

216 

HS ULX 

NGC4214 XMMl 

0.26699 


21 

LS IMBH cand. 

NGC4258 XMM3 

0.46503 


36 

LS IMBH cand. 

NGC4395 XMMl 

9.04609 

2.94683 

23 

HS ULX 

NGC4449 XMM2 

2.48586 


191 

LS IMBH cand. 

NGC4490 XMMl 

16.8513 

3.21972 

25 

HS ULX 

NGC4490 XMM2 

7.36612 

4.51554 

35 

HS ULX 

NGC4490 XMM3 

240.653 

176.04 

1354 

HS ULX 

NGC4490 XMM4 

1.66829 


128 

LS IMBH cand. 

NGC4490 XMM5 

12.7136 


978 

LS ULX 

NGC4631 XMMl 

10.6527 

8.59661 

66 

HS ULX 

NGC4736 XMMl 

31.6561 

27.3664 

211 

HS ULX 

NGC4945 XMM3 

0.44985 


35 

LS IMBH cand. 

NGC5204 XMMl 

22.4756 

2.20492 

17 

HS ULX 

NGC5204 XMM2 

5.57769 


429 

LS ULX 

M51 XMM2 

4.25898 

3.57502 

28 

HS ULX 

M51 XMM3 

2.10133 


162 

LS IMBH cand. 

M51 XMM4 

2.56064 


197 

LS IMBH cand. 

M51 XMM6 

46.7642 

39.5189 

304 

HS ULX 

M101 XMMl 

8.04916 

7.68224 

59 

HS ULX 

MIDI XMM2 

7.54268 

4.96709 

38 

HS ULX 

M101 XMM3 

10.9792 

1.03659 

8 

HS ULX 

NGC5408 XMMl 

20.9211 

11.5369 

89 

HS ULX 

Circinus XMMl 

70.3579 

56.1033 

432 

HS ULX 

Circinus XMM2 

208.746 

0.69929 

5 

HS ULX 

Circinus XMM3 

771.157 

0.212957 

2 

HS ULX 


^estimation of the bolometric luminosity, determined with an exponential 
cut-off in the power law at high energy (see text) 

b bolometric luminosity estimate for high-state ULXs where the power law is 
cut at twice kT* n (see text); units for both luminosity measurements in 10 39 
erg s~ l 
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c mass computed for objects radiating at 0.1 X L^dd (low-state objects) or 
L Edd (high-state objects; using L^), in units of Mq 

d classification based on the criteria set forward in the text: high-state ULX 
(HS ULX) , low-state ULX (LS ULX) , and low-state IMBH candidate (low-state 
object with lO 38 ergs - ' 1 < L b 0 i < 3 x 10 39 ergs -1 ) 




Hubble Type 

Fig. 1. — Distribution of galaxies by Hubble type among our archival XMM-Newton sample of 
nearby (< 8Mpc) galaxies. Our sample consists solely of spirals and irregulars. 
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i°g(f„/f„pt) 


Fig. 2.— The distribution of f x /fo P t for the brightest possible optical point source within the 
XMM-Newton error circle. We define f x as the unabsorbed X-ray flux in the 0.3 — lOkeV range and 
f opt as the optical flux obtained from the U filter of XMAfs OM (as described in text). These ratios 
do not represent the actual f x /f 0 pt of the sources but are an estimate of the minimum possible 
value. A majority of the sources had no optical point source within the X-ray contour and thus 
have ratios of f x / f opt far higher than those indicated in the plot. 




a 

D 

«-i 

o 

0 ) 

£ 

S 

P 

z 


Spectral Index (F) 

Fig. 3. — Distribution of the spectral indices (F) for low-state (solid) and high-state (dashed) 
objects. For Galactic low-state objects, typically F « 2.0, similar to our sample, while the high- 
state objects have a steeper F (McClintock & Remillard 2004), 
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Spectral Index (F) 


Fig. 4. — Relationship of luminosity vs. spectral index for low-state (rectangle) and high-state 
(triangle) objects. Sources represented by a circle are those where the Ax 2 value between the 
two-component and power law fits was very small. As expected from observations of Galactic 
stellar-mass black hole systems (McClintock & Remillard 2004), the classified low-state ULXs in 
our sample have, on average, lower X-ray luminosities than the corresponding high-state ULXs. 
We plot the mean values for both high-state and low-state objects with errorbars indicating the 
root mean square deviation. The outlying objects with spectral indices greater than 3.5 were not 
included in the mean or deviation calculations. 
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log(kT) log(kT) 


Fig. 5. — (left)Distribution of the blackbody temperature for high-state objects. (right) Relationship 
of blackbody temperature vs. luminosity (in the 0.3-10 keV band) for high-state objects. We see 
two peaks arise in the distribution, one centered around kT?s 0.1 and another at kT« 1. The peak 
with a low disk temperature also corresponds to the highest luminosities, suggesting that these may 
be high-state IMBHs. The sources with higher disk temperature also have lower luminosities. The 
spectra of these sources were also well fit by an inverse Comptonization model (a model succesfully 
used to fit some of the Galactic black hole X-ray binaries in the very high state). 
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Fig. 6. — Color-Color Diagram plotting soft vs. hard colors, as outlined in Done & Gierlinski 
(2003), for low-state (rectangle) and high-state (triangle) ULXs. The filled circles represent the 
sources with low A% 2 values between the two-component and power law models. A large number 
of our sources lie in the same range of this graph as the black hole sources examined by Done & 
Gierlinski (2003) (near the power law distribution, indicated by the solid line). The dashed line 
represents the color-color plot for a multi-colored disk model with different disk temperatures. The 
sources approaching this line were those well-fit by the Cornptonization model. Done & Gierlinski 
(2003) had no black hole sources in this region, but atolls and Z-sources, which were also well-fit 
by Cornptonization models. 
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log(4l 8 ) log(4lH) 


Fig. 7.— (left) Relationship of the far-infrared luminosity, as an indicator of star formation rate, vs. 
the number of ULXs for each galaxy. If ULXs are associated with star formation, we naively expect 
that the higher the FIR luminosity the more ULXs the galaxy will host, (right) The distribution of 
average number of ULXs / L fir bin for spirals (solid line) follows this expectation. The distribution 
of irregulars (dashed line) is not so easily interpreted. The numbers at the top indicate the number 
of spirals/irregulars in each of the luminosity bins. More irregulars would need to be included in 
this survey for meaningful statistics on this group. 






log(n H ) 

Fig. 8.— Distribution of the hydrogen column densities of ULX sources. The n# values were 
obtained through spectral fits using the wabs model in XSPEC. Galactic column densities towards 
the host galaxy were subtracted from the spectral fit values. A majority of our ULX sources have 
high column densities (> 10 21 cm -2 ), suggesting that some of this absorption originates with the 
local ULX environment. Bins to the left of the dashed line represent sources with column densities 
very close to the Galactic value and thus a simple subtraction is not statistically representative of 
the true value. 
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A. Spectral Simulations 
A . 1 . Two-component Model 

In order to determine the number of counts required to distinguish whether a blackbody 
component is statistically significant for the sources fit with a two-component model, we simulated 
spectra based on that of some of the brightest sources. We chose to simulate spectra of bright two- 
component spectra exhibiting three different cases: (I) the flux from the blackbody dominates over 
the power law component at 2 kT, (2) an intermediary case, and (3) the flux from the power law 
dominates over the blackbody component at 2 kT . Such simulations would allow us to determine the 
uncertainty in our claims of a combined fit being a better descriptor of the data. This is necessary 
because there is no a priori model which predicts the relative intensities of the two components and, 
as we know from studies of Galactic black holes, these components show a wide variety of relative 
intensities. To this end, we simulated spectra using the best-fit absorbed blackbody and power 
law model with the fakeit command in XSPEC. We chose (1) NGC 247 XMM1, (2) NGC 5408 
XMM1, and (3) Holmberg II XMM1 as our seed observations. These objects all have very high 
signal-to-noise and thus the fits are robust. The respective ratios of powerlaw flux to blackbody flux 
contributions at 2kT are: (1) << 1.0, (2) 1.77, and (3) 3.52. All of these sources have comparable 
blackbody temperatures indicative of our high-state ULX candidates (roughly kT ~ 0. 15). 

We simulated 500 spectra each, using the two-component model, for each of 2000 counts, 
1000 counts, 400 counts, and 200 counts for the PN. Each simulated spectrum, based on the best-fit 
blackbody and power law model, was fit with an absorbed blackbody and powerlaw model as well 
as an absorbed pure-powerlaw model. We placed the constraint that the blackbody temperature 
must lie within the range of 0.07 — 4.0 keV (the range at which it would be detectable in the XMM- 
Newton bandpass). We allowed the power law index to vary over the range 0 — 4 for the power 
law component of the combined blackbody and power law model. However, we placed a constraint 
that the power law component must lie within the range F .== 1.5 — 2.0 for the simple power law 
model to be consistent with our fits to the sources we claim are best fit by simple power laws. This 
constraint ensures that the spectral index would exhibit that of our classified “low-state” objects. 

When analysing and classifying our real spectra, we declared a detection of the thermal disk 
component if the addition of this component (to a baseline powerlaw model) led to an improvement 
of the goodness of fit parameter by at least A% 2 = 8. We can use the above simulations to address 
the detectability of a thermal disk component using this Ax 2 threshold as a function of the relative 
strength of the thermal component and the number of counts in the spectrum. For each simulation, 
we fit the spectrum with both a single absorbed power-law and a 2-component powerlaw and thermal 
disk model and compute the quantity Ay 2 = Xpow — Xp 0W +disk- Fig. Al, we plot the distribution 
of Ax 2 from our 500 simulations for the weak and strong blackbody component for spectra with 
400 counts and 2000 counts. It is clear that we cannot detect a weak thermal component in a 400 
count spectrum — the vast majority of the simulations 82%) result in Ax 2 < 8. However, even 
a weak thermal component is easily detected in a 2000 count spectrum (not a single simulation 
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gave Ax 2 < 8). The strong black body case is detectable with high significance even in a 400 count 
spectrum (< 1% of the simulations resulted in A% 2 > 8). 

When we increase the upper limit of the range of the spectral index in XSPEC for the power law 
model to V — 3.0, our confidence levels decrease. For a weak thermal component with 400 counts, 
all of the simulations result in Ax 2 < 8. At 2000 counts, only ~ 30% of the simulations for a weak 
thermal component have Ax 2 < 8. For a strong thermal component, 20% of the simulations yield 
Ax 2 < 8 for 400 counts while, as was the case for F == 2.0, none of the simulations gave Ax 2 < 8 
for 2000 counts. 

Thus allowing the F to “float freely” or remain unconstrained further decreases the confidence 
levels. When the F parameter is allowed to float the spectra are fit with higher T values in order 
to compensate for the missing blackbody component. When the upper limit was instituted at 2.0 
or 3.0, we found that all 500 simulated spectra were fit with a T = 2.0 or 3.0, respectively, for the 
simple power law model. This same affect is not seen in the T of the two-component model, where 
the value ranges between 1.4 and 4.0 with a peak in the distribution near that of the original model 
used to simulate the spectra. Thus, higher power law indices (> 3.0) can indicate the necessity 
of an added blackbody component. Since there are no Galactic black holes whose broad band 
-spectra are well fit by steep power laws it seems that restricting the allowed power law indices 
is more consistent with the nature of Galactic black holes. In fitting our sources with the three 
“standard” models, we allowed the F to float, thus the problem of a missing blackbody component 
being compensated for by a steep power law should not have factored into our classification criteria. 
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Ax 2 value 


Fig. A9.— Distribution of A% 2 from simulations. The Ay 2 values represent the difference between 
the unreduced x 2 of the absorbed power law model and the combined blackbody and power law 
model. The left panel shows the results for the “weak” blackbody component relative to power law 
for 500 simulated spectra at 400 counts (top) and 2000 counts (bottom). The right panel shows the- 
results for the “strong” blackbody component. For a weak blackbody component at 400 counts, 
the thermal component is undetectable. However, it is able to be distinguished as the number of 
counts is increased. A strong blackbody component is easily distinguishable at 400 counts. 
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A. 2. Simple Power Law Model 

Our next set of simulations sought to determine our confidence in the simple power law fit 
being an adequate descriptor of the spectra. Binning provides a problem in distinguishing between a 
powerlaw and a curvature in the spectrum at the low energy range, since the binning procedure can 
wash out a low kT blackbody from the spectrum. For this reason, we chose to simulate unbinned PN 
spectra for a source we categorized as a low-state object, IC 0342 XMM1 . IC 0342 XMM1 represents 
characteristics typical of our low-state candidates, namely it is within the proper luminosity range 
and it has a power law index of T & 1.7 (the median of the distribution for low-state objects is 2.03) 
and a hydrogen column density near the median of all the fitted values (where the median value 
is « 3 x 10 21 cm -2 and the value of IC 0342 X-l’s column density is 5.8 x lO 2 * cm -2 ). We chose 
this source for these reasons and the high number of counts in its PN spectrum. Instead of using 
the x 2 statistic (used for binned data), we chose to use the maximum-likelihood statistic, C-stat , 
in XSPEC (which uses unbinned data). 

We simulated 500 spectra using the best-fit parameters for the simple power law fit using 
the fakeit command. We fit the simulated spectra with two models: the simple powerlaw and a 
combined blackbody and power law model and computed the change in the goodness of fit, Ay 2 . 
For the two-component model, once again we placed the constraint that the blackbody temperature 
remain in the range that it would be detectable by XMM-Newton , 0.07 - 4keV. We followed this 
procedure for 1000, 2000, and 4000 counts. At both the 1000 and 2000 count level, the addition 
of a second component has no effect on the C-statistic (the distributions in C-space for both the 
power law fit and the two-component fit are indistinguishable). 

At the 4000 count level, the C-statistic distributions for the two models separate such that 
there is an 26% confidence that the two-component model is a better fit to the data. We find, 
when we examine the model parameters, that the power law index (F) for the two-component 
model ranges between 1.19 — 2.12 with the mean value « 1.74. The mean value for all three count 
levels used clustered around this value, though the range in F increased as the counts decreased. 
The mean blackbody temperature for simulations with 4000 counts was 1.17keV (with the range 
varying between the amount previously noted) with a median of 1.06 keV. For simulations with 
lower counts, the blackbody temperature becomes higher (1.24 keV for 1000 counts) with a higher 
median (1.83 keV for 1000 counts). This tells us that the fitting procedure tends to approximate a 
pure power law spectrum as a two-component spectrum with T equal to that of the true spectral 
index but with a blackbody temperature higher than those observed in our study (1.1 keV or higher) 
which in the XMM band can be approximated as a power law. If we found spectra in our sample 
that were best fit with a low spectral index and a high blackbody temperature, we might suspect 
that the spectrum’s true nature is a power law. We also note that if the hydrogen column density 
is large, much > 3 x 10 21 cm -2 , a low temperature (kT) blackbody component can be much more 
difficult to detect. 
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B. Additional Spectral Fits 

The following sources were not best fit by the standard models employed in this study: 


B.l. NGC300 XMM4 

This source was classified as a super-soft X-ray source by Kong & DiStefano (2003). We find 
that the standard single-component absorbed blackbody model is a much better model for this 
spectrum. In fact, the power law, bremsstrahlung, and combined models do not fit the data within 
the 90% confidence range. Fitting an absorbed blackbody, we find the best fit corresponds to the 
following parameters: njj = 1.38^55 x 10 21 cm~ 2 , kT = 0 . 059 t°;oo 5 k e V> anc * X 2 /d°f — 74.5/45. 
This fit yields an unabsorbed flux of 3.3 X 10~ 13 ergcm~ 2 s -1 . 


B.2. NGC4631 XMM4 

The spectrum of this source clearly identifies it as a super-soft X-ray source. As with NGG300 
XMM4, the standard models employed in this study did not adequately match the data. The best 
fitting model corresponds again to an absorbed blackbody. The corresponding parameters are as 
follows: %h = 6.2l°;g 6 x 10 21 cm~ 2 , kT = O.OTIq;^ keV, and X 2 /dof — 142.3/74. This fit yields an 
unabsorbed flux of 9.5 x 10~ 12 ergcm -2 s -1 . The position of this source shows it to be coincident 
with a globular cluster associated with that galaxy. This source was identified as a bulge X-ray 
source, possibly powered by accretion, in a ROSAT study of NGC4631 (Vogler k, Pietsch 1996). 


B.3. NGC4631 XMM5 

The spectrum of this source was best fit with an absorbed power law -f an absorbed vapec 
model. This indicates the prescence of hot gas, indicating a possible thermal X-ray source. 


B.4. NGC4945 XMM5 

The spectrum of this source was not adequately fit with any of the standard models used in 
this investigation. The spectrum exhibits a prominent Fe K line in the PN spectrum that is well fit 
by a gaussian ( zgauss ) at 6.4 keV. We find that the entire spectrum is best fit with a partial covering 
fraction absorption model ( pcfabs ) in combination with the normal absorption, a power law, and 
a gaussian. The best fit parameters yield: absorption column density, njj = 1.79 X 10 2! cm -2 , 
partial covering absorption, nn ~ 18.4 x 10 21 cm” 2 , partial covering fraction = 0.82, P = 1.6, and 
X 2 /dof = 61.8/57. The source is clearly located within the optical galaxy, and is thus unlikely to 
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be a background AGN. 


B.5. M51 XMM5 

The spectrum and luminosity (L.y « 1.9 X 10 42 ergs -1 ) of this source suggests that it is an 
AGN. The location of the source, from the Digital Sky Survey, places it within the dwarf companion 
of M51 making a value of the optical flux hard to constrain. The best fit to this source was an 
absorbed blackbody + power law and the spectral parameters are listed in Table 4. 


B.6. M83 XMM2 

Like NGC4945 XMM5, this source was best fit by a partial absorption model. However, this 
source showed no evidence of an Fe K line. We fit this source’s spectra using a partial covering 
fraction absorption model in combination with the normal absorption model and a power law. 
The best fit parameters yield: absorption column density, n# = 2.1 x 10 21 cm -2 , partial covering 
absorption, n# = 43.5 x 10 21 cm -2 , partial covering fraction =s 0.86, T = 2.95, and X 2 / dof = 
83.5/84. The unabsorbed flux in the range of 0.3-10 keV equals 1.37 x 10 -12 ergcm -2 s -1 . 


B.7. Inverse Compton Scattering Sources 

Table C4 includes the parameters for the “ULX” sources best fit by the compST model. A 
discussion of these sources and interpretation of the data is included in section 4.3. 


C. Additional Tables 
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15 Unable to model sprectrum due to an error in x 2 -space 

16 bright sources are too close, but resolvable by Chandra (Summers et al. 2004) 
17 alternate ID from Jenkins et al. (2004) 



bright sources near the nucleus are unresol vable, but seen by Chandra (Smith & Wilson 2001) 
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Table C7. Bright, Identifiable Background and Foreground Sources 


Galaxy 

HA (h m s) 

Dec (o / //) 

Identification 

NGC 247 

0 46 51.7 

-20 43 30 

QSO B044-2059 

NGC 300 

0 55 26.7 

-37 31 25.6 

HD 5403 (Star) 

NGC 625 

01 34 42,4 

-41 36 15.2 

QSO B0132-4151 

NGC 1569 a 

04 31 16.9 

+64 49 50 

CXOU J0431 16 .8+644950 (Star) 

NGC 1569 

04 31 14.2 

+64 51 07.9 

CXOU 043114.0+645107 (Star) 

NGC 1569 

04 31 25.4 

+64 51 53.8 

CXOU 043125.1+645154 (AGN) 

NGC 1705 

04 54 01.2 

-53 21 12.3 

WGA J0454.0-5320 (M star or elliptical galaxy) 

NGC 2403 

07 35 09 

+65 40 27.5 

HD 59581 (Star) 

NGC 4258 

12 IS 08.9 

+47 16 08.3 

QSO J1218+472 

M83 

13 36 45.6 

-29 59 13.9 

2MASX J 13364579-2959122 (Galaxy) 

M83 

13 36 13.9 

-29 56 13 

RX J 1336 15- 2957.8 (Galaxy) 

NGC 5253 

13 39 50.6 

-31 34 11.1 

CD-30 10790 (Star) 

M101 

14 02 30 

+54 21 18.2 

[WIP99] H13 (Star) b 

NGC 5408 

14 03 27.5 

-41 25 18.5 

(Star) 


identification for objects in NGC 1569 from Martin, Kobulnicky, & Heckman (2002) 
b confirmed by K. Kuntz using HST ACS 



Table C8. XMM-Newton power law fit for best fit two-component spectra 


Source n# a P x 2 /dot Fx c 


112.2/95 1900 2200 


NGC247 XMMl 
NGC253 XMM1 

NGC253 XMM2 (obs 2) 
NGC253 XMM3 

NGC253 XMM4 

NGC253 XMM5 

NGC253 XMM6 
NGC253 XMM7 
NGC300 XMM1 
NGC300 XMM2 
NGC300 XMM3 
NGC300 XMM5 
NGC300 XMM6 
NGC1313 XMM1 
NGC1313 XMM2 
NGC1313 XMM3 
ICO 34 2 XMM3 
NGC1705 XMMl 
NGC1705 XMM2 
NGC1705 XMM3 
NGC2403 XMMl 
NGC2403 XMM2 
NGC2403 XMM3 
Holmll XMMl (obs 1) 
Holm I XMMl 
MSI XMMl 

M81 XMM2 
M81 XMM3 
M81 XMM4 
M81 XMM5 
Holm IX XMMl 
Sextans A XMMl 
NGC4214 XMM2 
NGC4258 XMMl 
NGC4258 XMM2 (obs 1) 
NGC4395 XMMl 
NGC4395 XMM3 
NGC4449 XMMl 
NGC4449 XMM3 
NGC4490 XMMl 
NGC4490 XMM2 
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i-»7±S:iS 

1 c+0-07 
_ 0 07 

2-61+2:21 


9 04+°* i4 
0 o 7 

o 9 + 0.07 
■ "^—0.07 

9 aq+ 0-02 

4.uy_o.o 2 

o n+o.3 
^ ■ —0.3 

i 7Q+0-07 
A * ry — 0.07 

7.3 

6.13 

0 97 + 0 * 25 
u * y/ -o.4i 

■1.58+2:12 

?+? 
■ _? 

0 . 88 + 2:11 

i-oioi 

i.53+2:ii 

t y+0.08 
A *'* —0.08 

1 84+ 0 * 03 
i,54 -0.03 

O-l' 8 +o-.S 

9 oc + 0.12 

— 0.07 

0,3+23 

9 ao+0.43 


1 q+ 0.14 
A * y — 0.13 

o e+0.9 

1 QQ+0.16 

.0.15 

o 7+0.5 

4.93+232 

? + ? 

* — ? 

148+2:22 

g o+ 0 . 9 

9 99+O.I4 
Z * ZZ -0.12 

o o+D.5 
*•*-0.4 

3 3g + 0 - 29 
— 0.23 

0 g3 + 0 - 14 

U. 50 — 042 

9 co+0.17 
^• 00—0 16 

6 3 + 1 ' 3 
■ *— 1.0 

2 3g +°* 19 
■ .r-Q.,17 


262.6/232 

3.1 

3.6 

611.6/582 

2.9 

3.3 

507.4/500 

1.2 

1,4 

91.9/83 

0.73 

1.2 

381.8/409 

0.98 

1.6 

73.6/59 

0.52 

0.85 

321.4/293 

0.29 

.0,48 

31.8/25 

0.32 

0.53 

283.8/298 ‘ 

1.1 

1.3 

435/409 

0.93 

1.5 

357/342 

1.3 

2.2 

469.8/422 

0.81 

0.6 

133.98/99 

0.27 

0.21 

101.9/81 

0.19 

0.15 

54.2/56 

0.14 

0.10 

47.6/37 

0,06 

0,15 

219.8/203 

0.42 

0.88 

464.1/421 

1,9 

4.0 

778.3/426 

3.3 

6.9 

185,8/109 

1.7 

3.1 

61.9/88 

0.12 

0.37 

91/76 

0.078 

0.24 

80.9/67 

0.17 

0.53 

92.2/81 

2.2 

3.6 

179.5/151 

1.3 

2.0 

92.6/107 

0.81 

1.3 

1134.2/976 

12 

10 

102.8/95 

0.48 

1,7 

1849.9/1245 

4.5 

7.0 

224.9/208 

4.3 

6.7 

1358.2/618 

48.5 

75,2 

81.35/80 

2.5 

3.9 

66.4/52 

0.35 

0.54 

97.5/82 

0.44 

0.68 

1000.9/882 

9,4 

15 

271.4/275 

0.56 

0.13 

50.9/46 

0.16 

0.14 

101.4/78 

0.06 

0.04 

97.5/63 

0.43 

2.6 

195.1/156 

7.2 

14 

55.9/58 

0.25 

0.48 

103/118 

1.2 

1.36 

154/89 

1.3 

1.5 

101.5/65 

1.2 

8,7 

49.5/56 

0.92 

6.7 
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Table C 8 — Continued 


Source 


r 

X 2 /dof 

Fx c 

Lx d 

NGC4490 XMM3 

9 4 + 1 ' 5 
y * 4 - 1.2 

9 qc: + 0-24 
z * yD ~ 0.20 

76.7/80 

1.5 

11 

NGC4631 XMMl 

9 o + 0*16 

2 *°-0.15 

9 1 0 +O.O 6 

383.4/347 

0.76 

5.1 

NGC4631 XMM 2 

1 q+°- 4 
i ' y ~0.3 

AUi ~ 0.12 

119.5/99 

0.23 

1.5 

NGC4631 XMM3 

0.63^1 

i co + 0.1 
08 

146/98 

0.15 

1.0 

NGC4736 XMMl 

0.95±g;| 

2 02 + °* 26 
2.U2„o.25 

62.8/53 

0.36 

0.80 

NGC4945 XMMl 


1 oc+0-08 

J-*OO _ 0 jo 

116/122 

0,9 

1.0 

NGC4945 XMM 2 

3 4 +°- 6 

i co+0.09 

J..OO_Q iq 

114.5/115 

0.71 

0.82 

NGC4945 XMM4 

5-2+S't 

9 cq+0 - 19 

2 *oy_ 0 .i7 

75.5/62 

0.49 

0.56 

NGC5204 XMMl 

o-siiol 

9 11 +0*0 , 4 
■ ■ — 0.04 

592.1/561 

2.0 

5.5 


1 1+ 01 
i ' i — 0 . 1 ' 

9 A-t +0.07 
■ ■ — 0.07 

533/498 

3.0 

8.3 

M51 XMMl 

, ,40.30 
i - 1 ~0.27 

9 67+0*' 20 

Z.O/ _ 0 .16 

110.5/82 

0.34 

2.8 

M51 XMM 2 

9 o+O.SO 
z *°-0.30 

2 50 ^"0-22 

75.2/70 

0.52 

3.3 

M51 XMM5 

2.7 

3.08 

256.0/72 

0.43 

2.7 

M51 XMM 6 

0 n+0,83 
Z - U ^0.72 

2 50"*"°' 33 
z.su _ 0 25 

40.97/43 

0.13 

0.83 

M51 XMM7 

n c+0.39 
-0.46 

i qcj + 0.23 

37.8/31 

0.11 

0.66 

M83 XMMl 

1 Q +°‘ 34 
• ■ — 0.31 

9 09+6.13 
**^-0.12 

210.9/211 

0.64 

2.8 

M83 XMM4 

e-oil;! 

2 54+O -26 
■ ■ — 0.23 

91.4/91 

0.4 

1.8 

MIDI XMMl 

0.56l»;l! 

1 QR+ 993 

l* y °_ 0 .08 

303/233 

0,45 

2.9 

M 101 XMM 2 

9 9+0.25 
Z — 0,23 

1 R5 + 007 

288.8/263 

0.81 

5.3 

Ml 01 XMM3 

i 5+0 .4 
i * £> -0.3 

2 70+° ' 2i 
Z * /U -0.17 

148.9/133 

0.51 

3.4 

M 101 XMM4 

9 9+0.45 

tS- 42 

2 - 25 to.i° 

165.7/140 

0.38 

2.5 

M 101 XMM5 

1 O+ 0.2 

— 0.3 

2.28l®;l! 

0 57+0.12 

47.9/46 

0.12 

0.8 

NGC5408 XMMl 

1 6 + 0-2 

396.8/339 

7.04 

19.4 

CIRCINUS XMMl 

7 - 6 lS'i 

9 1 5+0.05 
Z * i0 -0.04 

762.9/863 

4.6 

8,8 

CIRCINUS XMM 2 

n.7l!;! 

0.06 

517.9/432 

2.7 

5.2 

CIRCINUS XMM3 

9.011;! 

9 57+O.4O 
Z * &/ -0.17 

285.2/262 

0.32 

0.61 


a total column density in units of 10 21 cm” 2 

^improvement in x 2 over the single-component power law model 

c unabsorbed flux in the 0.3-10 keV band in units of 10“ 12 erg cm” 2 s“* 

d unabsorbed luminosity in the 0.3-10 keV band, using the distances 
quoted in Table 1, in units of 10 39 erg s” 1 

e absorption level frozen at Galactic level 
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Table C9. Best-Fit Absorbed Comptonization Model Parameters 


ID 

nH a 

kT b 

tau c 

x 2 

Fx d 

NGC 253 XMM2 

1 Q + 0.08 

i-2S±g;J| 

19.59+16 

464/498 

1.47 

NGC 2403 XMMl 

i-95±J;l 

0.981°;}® 

25.4l7;| 

82.8/85 

1.4 

NGC 4490 XMMl 

4-7±|;l 

o-96±°:il 

27.0±J 8 9 3 

66.5/64 

0.66 

NGC 4490 XMM2 

5-0±il 

■t 9*1+0,21 

18.8.11;® 

45.7/55 

0.67 

MIDI XMM2 

1 £> + 0.26 
- 1 *t>-0.24 

i 94+0-09 

23.3l®j 

256/262 

0.65 

M101 XMM3 

, 1 +0.43 
i *1-0.40 

1.13 

15.2 

128/132 

0.41 

Circinus XMM2 

8.8+S1 

0.821855 

29.7l 8 J 

437.2/430 

0,5 

Circinus XMM3 

6-7 tli 

0.93l°-|® 

28.1,1^ 

273.1/261 

0.17 


a total column density in units of 10 21 cm -2 
b temperature in keV 
c optical depth 

d unabsorbed flux in the 0.3-10 keV band in units of 10” i2 erg cm -2 s _1 
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Table CIO. XMM-Newton Galaxy Observations 


Galaxy 

S6o (Jy) 

Sioo (Jy) 


LFJH b 

No. of ULX 

NGC247 

7.93 

27.32 

0.602 

0.687 

1 

NGC253 

998.73 

1861.67 

55.92 

93.10 

3 

NGC300 

23.08 

74.45 

1.688 

1.324 

0 

NGC625 

5.09 

9.08 

0.280 

0.230 

0 

NGC1313 

35.97 

92.00 

2.329 

4.845 

2 

IC0342 

255.96 

661.68 

16.66 

30.32 

3 

NGC1569 

45.41 

47.29 

2.072 

0.635 

0 

NGC1705 

0.970 

2.580 

0.064 

0.199 

0 

MRK 71 

3.51 

4.67 

0.173 

0.239 

1 

NGC2403 

51.55 

148,49 

3.547 

5.378 

2 

Holmberg II 

1.15 

2.62 

0.070 

0.061 

1 

Holmberg I 





1 

M81 

44.73 

174.02 

3.647 

5.655 

1 

M82 

1271.32 

1351.09 

58.35 

106.2 

1 

Holmberg IX 





1 

Sextans A 

0.255 

0.674 

0.017 

0.004 

0 

IC 2574 

2.41 

10.62 

0.212 

0.329 

0 

NGC 4214 

17.87 

29.04 

0.947 

0.826 

1 

NGC 4258 

21.60 

78.39 

1.690 

10.48- 

1 

NGC 439 5 

4.21 

12.90 

0.299 

0.573 

1 

NGC4449 

37.00 

58.28 

1.937 

2.199 

1 

NGC4490 

47.79 

85.94 

2.636 

19.19 

5 

NGC4631 

82.90 

208.66 

5.324 

35.83 

1 

NG C4736 

62.41 

135.34 

3.734 

8.261 

4 

NGC4945 

588.11 

1415.5 

36.95 

42.49 

0 

NGC 5204 

2.33 

5.35 

0.143 

0.395 

2 

M51 

108.68 

292.08 

7.213 

44.74 

5 

M83 

266.03 

638.63 

16.69 

76.79 

0 

NGC 5 253 

30.00 

30.92 

1.365 

16.72 

0 

Ml 01 

88.04 

252.84 

6.048 

39.63 

4 

NGC 5408 

2.825 

2.958 

0.129 

0.356 

1 

Circxnus 

248.7 

315.85 

12.06 

23.10 

4 


a flux in units of 10 9 erg cm 2 s -1 
b far- infrared luminosity in units of 10 42 ergs -1 * 



